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1.  EXECUTIVE  SUMMARY 

Cryogenic  microcoolers  was  characterized  for  optimum  thermal  performance  and 
thermo-mechanical  reliability  when  employed  for  high  power  laser  diodes  (LDs).  The  goal  was 
achieved  by  conducting  two  specific  tasks:  Development  of  Cryogenic  Microcooler  and 
Apparatus  for  High-Power  LD  Bar  (Part  I)  and  PoF-based  Reliability  Assessment  of 
LD/Microcooler  Subassembly  Subjected  to  Cryogenic  Operating  Conditions  (Part  II). 

In  Part  I,  an  open-loop  liquid  nitrogen  cooler  for  the  cryogenic  operation  of  LD  bars  was 
developed  and  tested.  During  design  of  the  cooler,  consideration  was  paid  to  the  needs  and 
conditions  under  which  a  field-deployed  system  would  operate.  While  an  open-loop  design 
presented  particular  challenges  for  the  precise  characterization  of  cooler  performance,  it 
dispensed  with  the  bulk,  complexity,  and  power  consumption  associated  with  a  closed  nitrogen 
refrigeration  system.  In  keeping  with  this  design  direction,  modifications  to  the  test  loop  were 
made  to  reduce  system  startup  time  and  prevent  ice  buildup  on  the  LD  emitter  surface.  These 
included  an  aerogel  insulation  layer  and  nitrogen  shield  gas  sourced  from  the  LN2  reservoir,  as 
opposed  to  isolating  the  test  loop  within  a  high- vacuum  chamber.  An  internally-routed  bypass 
system  demonstrated  a  reduction  in  the  time  for  cool-down  from  a  room  temperature  to  as  little 
as  four  minutes,  compared  to  previous  hour-long  startups.  These  resulted  in  compact,  coolant- 
efficient,  and  responsive  system.  The  results  from  cryogenic  testing  of  a  mounted  LD  bar  showed 
that  the  cryogenic  cooler  and  its  housing  added  only  20%  additional  thermal  resistance  to  that  of 
the  commercial  block-mounted  LD  bar.  The  base  of  this  block  mount  was  brought  to  -180°C 
under  a  dissipated  power  load  of  55  W. 

In  Part  II,  the  reliability  of  LD/microcooler  subassemblies  subjected  to  cryogenic  operating 
conditions  was  assessed  by  proposing  a  hierarchical  PoF-based  reliability  assessment  model.  In 
the  model,  three  aspects  of  the  reliability  were  considered:  optical  reliability,  LD  bar  reliability, 
and  thermal  reliability.  The  junction  temperature  distribution  measurement  was  essential  for  the 
model.  A  hybrid  experimental/numerical  method  was  proposed  and  implemented  to  predict  the 
junction  temperature  distribution  of  a  high  power  LD  bar  first.  Then,  a  novel  method  was 
developed  to  predict  the  spectral  power  distributions  (SPDs)  of  individual  emitters  by 
deconvoluting  the  SPD  of  an  LD  array.  This  method  was  utilized  to  predict  the  SPD  change 
caused  by  the  die  attach  crack  of  the  LD  bar.  Thermal  and  mechanical  models  were  subsequently 
developed  to  predict  thermo-mechanical  stresses  of  the  subassembly  and  thus  to  be  able  to 
predict  the  potential  failure  locations  and  the  crack  propagation  within  the  die  attach  during  test 
conditions.  Finally,  the  effect  of  die  attach  crack  on  the  SPD  change  was  calculated  for  the 
physics-of-failure  (PoF)-based  reliability  assessment. 
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2.  PART  I:  DEVELOPMENT  OF  CRYOGENIC  MICROCOOLER  AND 
APPARATUS  FOR  HIGH-POWER  LD  BAR 

Detailed  heat  transfer  data  for  two-phase  nitrogen  flow  in  complex  structures  (e.g.  microfluidic 
coolers)  is  underrepresented  in  literature.  With  little  prior  work  on  which  to  model  the  cryogenic 
cooler,  considerable  attention  was  paid  to  predicting  cooler  performance,  choosing  a  candidate 
design,  and  characterizing  the  fabricated  cooler  prior  to  integrating  the  LD  bar  for  final  testing. 
This  chapter  summarizes  those  efforts,  and  concludes  by  detailing  the  results  of  cryogenic 
operation  of  the  LD  bar. 

2.1  Development  of  Cryogenic  Microcooler 

2.1.1  Prediction  of  Pin-Fin  Thermo-fluid  Characteristics 

The  use  of  microcoolers,  relying  on  evaporation  and  two-phase  flow  through  arrays  of  micro  pin 
fins,  is  not  yet  established  in  the  electronic  cooling  community  and  little  research  has  been  done 
to  quantify  and  correlate  the  heat  transfer  rates  and  pressure  drops  for  such  configurations.  The 
available  correlations  include  no  LN2  data,  and  -  like  many  other  two-phase  flow  correlations  - 
rarely  provide  better  than  an  accuracy  of  +/-  30%.  Nonetheless,  these  correlations  can  still 
identify  general  thermofluid  trends  for  two-phase  flow.  Among  these,  the  pin-  fin  correlation 
derived  by  Peles  &  Krishnamurthy  [1]  has  been  extensively  evaluated  by  other  investigators  and 
found  to  provide  the  most  consistent  results.  Accordingly,  the  purpose  of  this  preliminary  study 
is  to  observe  general  parametric  trends  for  two-phase  flow  through  micro  pin-  fin  structures  with 
LN2  using  the  correlation  originally  derived  by  Peles  &  Krishnamurthy  [1]  and  modified  by 
Reeser  et  al.  [2]  and  subsequently  tuning  the  correlation  to  the  data  sets  obtained  in  the  current 
experimental  effort.  Of  specific  interest,  during  both  validation  and  actual  Laser  Diode  bar 
testing,  are  the  effects  of  parameters  such  as  flow  quality,  heat  flux,  mass  flux,  pin-fin  geometry, 
and  flow  regime  on  the  thermal  performance  (heat  transfer  coefficient)  of  the  final  micro  pin-fin 
cooler  design.  The  Reeser  et  al.  [2]  correlation  was,  thus,  used  to  gather  some  insight  on  these 
parametric  effects  and  to  guide  the  initial  design  of  the  micro  pin  fin  cooler. 
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The  Reeser  et  al.  [2]  correlation  is  given  in  full  detail: 


htp  —  C(<A  )  hSp 

where, 
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D/=<4,  =  hydraulic  diameter  of  single  pin 
SL  =  streamwise  pin  pitch 
Sr  =  traverse  pin  pitch 
N  =  total  number  of  pins 
hf=  height  of  pin 

dh  =  gap  above  pin  (0  for  this  study) 

Re  =  Reynolds  number 
Pr  =  Prandtl  number 

Prw  =  Prandtl  number  for  wall  temperature 

f  =  friction  factor 

kf  =  thermal  conductivity  of  fluid 

Xw  =  Martinelli  parameter 

x  =  flow  quality 

xe  =  flow  quality  at  exit 

G  =  mass  flux 

pf =  density  of  liquid 

pv  =  density  of  vapor 

hsp  =  single-phase  heat  transfer  coefficient 
htp  =  two-phase  heat  transfer  coefficient 
Cnu>  Cp  ^2’  ^3,  C4,  C5  =  empirical  constants 


As  shown,  the  correlation  takes  into  account  the  pin-fin  geometry,  fluid  &  thermal  properties, 
mass  flow  rate,  and  flow  quality.  The  range  of  conditions  tested  was  based  on  the  parametric 
space  in  which  the  correlation  was  derived.  The  baseline  conditions  are  always  plotted  in  black 
and  are  as  follows: 


•  LN2 

•  Tsat  =  80K 

•  800  kg/m2-s  mass  flux 

•  12x12  mm2  pin-fin  channel 

•  Inline  configuration  of  pins 

•  150  pm  wide  pins 

•  300  pm  pitch 

•  300  pm  height 

•  -275%  Area  Enhancement 

2.1.1. 1  Mass  Flux 

The  effect  of  mass  flux  on  the  heat  transfer  coefficient  vs.  flow  quality  profile  is  demonstrated  in 
Figure  1  for  the  baseline  conditions.  Note  that  the  base  heat  transfer  coefficient  takes  into 
account  the  surface  area  enhancement  from  the  pin-fin  structure.  As  shown,  the  correlation 
predicts  a  “softened”  M-shape  profile  and  a  rapid  deterioration  in  the  heat  transfer  coefficient 
(caused  by  dryout)  after  a  flow  quality  of  roughly  65%.  The  second  peak  in  the  profile  occurs  at 
flow  qualities  between  30  to  50%,  shifting  to  higher  qualities  with  decreasing  mass  flux.  This 
trend  is  attributed  to  the  higher  velocities  at  higher  mass  fluxes,  which  trigger  an  earlier 
breakdown  (dryout)  in  the  evaporating  liquid  film.  Overall,  increasing  the  mass  flux  enhances 
the  heat  transfer  coefficient  at  lower  qualities  and  decreases  the  heat  transfer  coefficient  at  higher 
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qualities,  with  a  transition  quality  that  decreases  with  increasing  mass  flux. 


Figure  1:  Influence  of  Mass  Flux  on  the  Base  Heat  Transfer  Coefficient  vs.  Flow  Quality 

Profile 

LNi,  Tsat  —  80K,  12  x  12  mm3  channel,  inline  configuration,  150  pm  wide  pins,  300  pm  pitch ,  and 

300  pm  height. 


2.1.12  Pin-Fin  Height 

The  effect  of  pin-fin  height  on  the  base  heat  transfer  coefficient  vs.  flow  quality  profile  is 
displayed  in  Figure  2  for  the  baseline  mass  flux  of  800  kg/m2-s.  As  shown,  the  heat  transfer 
coefficient  increases  across  the  entire  flow  quality  range  with  increasing  pin-fin  height.  The  local 
fin  wall  heat  transfer  coefficient  does  diminish  with  increasing  pin-fin  height  but  is  offset  by  the 
enhancement  in  surface  area  due  to  the  taller  fins,  resulting  in  an  overall  steady  increase  in  the 
base  heat  transfer  coefficient  with  fin  height.  According  to  the  correlation,  it  is  clearly 
advantageous  to  increase  the  pin-fin  height  above  the  300  pm  baseline,  thus  enhancing  the  heat 
transfer  coefficient  while  reducing  the  pressure  drop  as  a  result  of  the  increased  hydraulic 
diameter. 
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Figure  2:  Influence  of  Pin  Fin  Height  on  the  Base  Heat  Transfer  Coefficient  vs.  Flow 
Quality  Profile  for  a  Mass  Flux  of  800  kg/m2-s  and  Pin  Height  of  100, 200, 300, 400,  and 

500  pm 

LNi,  Tsat  =  80 K,  12  x  12  mm2  channel,  inline  configuration,  150  pm  wide  pins,  and  300  pm  pitch. 

2.1. 1.3  Streamwise  &  Traverse  Pin-Fin  Pitch 

The  effect  of  streamwise  and  traverse  pin-fin  pitch  on  the  base  heat  transfer  coefficient  vs.  flow 
quality  profile  is  demonstrated  in  Figure  3  and  Figure  4  for  the  baseline  mass  flux  of  800  kg/m2-s 
and  pin-fin  width  of  150  pm.  As  shown,  the  heat  transfer  coefficient  increases  across  the  entire 
flow  quality  range  with  decreasing  pin-fin  pitch.  Reducing  the  pin-fin  pitch  increases  the  wetted 
surface  area  (for  a  given  channel  footprint)  and  thins  the  liquid  film  on  the  fins  thus  enhancing 
the  base  heat  transfer  coefficient.  Additionally,  the  influence  of  streamwise  and  traverse  pin-fin 
pitch  on  the  heat  transfer  coefficient  is  identical  (i.e.,  Figure  4  and  5  are  identical). 

Accordingly,  it  is  thermally  advantageous  to  reduce  the  pin-fin  pitch,  but  thermal  gains  will 
come  at  a  cost  of  a  larger  pressure  drop  across  the  array. 
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Figure  3:  Influence  of  Streamwise  Pin  Fin  Pitch  on  the  Base  Heat  Transfer  Coefficient  vs. 
Flow  Quality  Profile  for  a  Mass  Flux  of  800  kg/m2-s 

LN2,  Tsat  =  80K,  12  x  12  mm2  channel,  inline  configuration,  150  pm  wide  pins,  and  300  pm 

height. 
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Figure  4:  Influence  of  Traverse  Pin  Fin  Pitch  on  the  Base  Heat  Transfer  Coefficient  vs. 
Flow  Quality  Profile  for  a  Mass  Flux  of  800  kg/m2-s 

LNi,  Tsat  =  80K,  12  x  12  mm2  channel,  inline  configuration,  150  pm  wide  pins,  and  300  pm 

height. 


2.1. 1.4  Variable  Pin-Fin  Width  with  Constant  Gap 

The  effect  of  pin-fin  width  -  for  a  constant  gap  of  150  pm  between  adjacent  pin-fins  -  on  the 
base  heat  transfer  coefficient  vs.  flow  quality  profile  is  demonstrated  in  Figure  5  for  the  baseline 
mass  flux  of  800  kg/m2-s  and  pin-fin  height  of  300  pm.  As  shown,  reducing  the  pin-fin  width 
significantly  increases  the  heat  transfer  coefficient  across  the  entire  flow  quality  range,  reflecting 
-  in  part  -  the  increase  in  wetted  surface  area  in  the  fin  array  and  -  in  part  -  the  thinner  film 
thickness.  However,  as  the  “minimum”  thickness  of  the  pin-fins  is  limited  by  machining 
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constraints,  (especially  for  high  aspect  ratio  pins),  the  baseline  pin  width  of  150  pm  was  chosen 
for  the  final  design. 


Flow  Quality,  X 


Figure  5:  Influence  of  Pin-Fin  Width  on  the  Base  Heat  Transfer  Coefficient  vs.  Flow 
Quality  Profile  for  a  Mass  Flux  of  800  kg/m2-s  and  Constant  Gap  of  150  pm  between 

Adjacent  Pins 

LNi,  Tsat  =  80K,  12  x  12  mm 2  channel,  inline  configuration,  and  300  pm  pin-fin  height. 

2.1.2  Summary  of  Preliminary  Correlation  Study 

According  to  the  results  obtained  by  exercising  the  Reeser  et  al  [2]  correlation  across  the 
parametric  space  of  interest  in  this  study,  it  is  thermally  advantageous  to  increase  the  mass  flux 
(for  low-  to  mid-range  qualities)  and  pin-fin  height,  and  decrease  streamwise/transverse  pin-fin 
pitch,  and  pin-fin  width.  These  are  the  main  trends  observed  from  this  parametric  study. 
However,  the  correlation  derived  by  Peles  &  Krishnamurthy  [1]  was  developed  using  relatively 
long  pin-fin  arrays,  while  -  in  the  present  application  -  the  pin-fin  array  is  significantly  shorter. 
Limited  chip-  scale  microgap  data  [3-4]  suggests  that  there  are  significant  differences  in  the 
thermofluid  behavior  of  short  and  long  microchannels.  It  is,  consequently,  important  to  obtain 
relevant  experimental  data  for  such  short,  pin  fin  channels  before  finalizing  the  design  of  the  LN2 
microcooler.  In  the  next  section,  data  from  such  an  experiment  using  a  fluid  similar  to 
LN2  (3M  FC-72)  are  compared  to  the  predictions  from  the  correlation. 

2.1.3  Experimental  Short-Array  Micro  Pin-Fin  Data 

In  the  previous  section,  a  two-phase  micro  pin-fin  correlation  derived  by  Reeser  et  al.  [1,  2]  was 
used  to  parametrically  analyze  the  thermal  performance  of  a  LN2  pin-fin  channel  for  various 
mass  fluxes,  flow  qualities,  and  pin-fin  geometries.  According  to  the  correlation,  it  is  thermally 
advantageous  to  increase  the  mass  flux,  pin-fin  height,  decrease  streamwise  and  transverse  pin- 
fin  pitch,  and  pin-fin  width.  However,  the  correlation  derived  by  Reeser  et  al.  was  developed 
using  relatively  long  pin-fin  arrays.  To  further  increase  confidence  in  the  selected  cooler  design 
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parameters,  the  correlation  predictions  were  benchmarked  against  a  short-array  pin-fin 
experimental  study  using  3M  FC-72  as  the  working  fluid. 

An  overview  of  the  pin-fin  channel  and  test  section  used  in  the  benchmarking  study  is  shown  in 
Figure  6  and  Figure  7,  respectively.  The  pin-fins  are  machined  out  of  the  copper  pedestal  shown 
in  Figure  6,  which  is  inserted  into  the  test  section  shown  in  Figure  7.  The  pin-fins  are  positioned 
in  an  inline  configuration,  are  150  pm  wide,  have  a  pin-to-pin  pitch  of  300  pm,  and  are  either 
100  pm  or  500  pm  tall.  The  entire  pin-fin  array  is  12  mm  x  12  mm  in  size,  amounting  to  a  total  of 
1600  pin-fins.  Heat  is  applied  to  the  bottom  of  the  pin-fins  with  a  ceramic  resistive  heater  and 
temperature  is  measured  using  five  embedded  thermocouples,  positioned  along  the  centerline  of 
the  pin-fin  channel  and  1  mm  below  the  surface.  The  pin-fin  channel  is  inserted  into  a  PEEK  test 
section,  which  channels  fluid  to  and  from  the  pin-fins.  The  upper  surface  of  the  pin-fins  is 
confined  with  a  sapphire  window,  which  is  used  for  top-down  photographic  and  mid-wave 
infrared  visualization. 


Figure  6:  Copper  Micro  Pin-Fin  Channel 
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Top  View 


The  experimentally  measured  base  average  heat  transfer  coefficient  and  pressure  drop  vs.  flow 
quality  profiles  are  shown  in  Figure  8  and  Figure  9,  respectively,  for  both  the  100  pm  and 
500  pm  tall  pin-fin  channels.  Data  for  the  100  pm  tall  pin-fin  channel  is  presented  at  three 
different  liquid  flow  rates:  0.5,  0.75,  and  1.0  mL/s  and  data  for  the  500  pm  tall  pin-fin  channel  is 
presented  for  one  flow  rate:  1.0  mL/s.  The  base  surface  area  enhancement  is  1.67x  and  4.33x  for 
the  100  pm  and  500  pm  tall  pin-fin  channels,  respectively.  The  predicted  base  heat  transfer 
coefficient,  using  the  correlation  derived  by  Reeser  et  al.  [1,  2],  is  plotted  as  a  dotted  line.  As 
shown,  the  heat  transfer  coefficient  and  pressure  drop  steadily  rises,  on  average,  with  increasing 
flow  quality,  heat  flux,  and  flow  rate.  The  100  pm  and  500  pm  tall  pin-fin  channels  reach  a 
maximum  heat  transfer  coefficient  of  55,000  and  43,000  W/m2-K  at  a  flow  quality  of  99  and 
94%,  respectively,  for  a  flow  rate  of  1.0  mL/s.  On  average,  the  base  heat  transfer  coefficient  for 
the  100  pm  and  500  pm  tall  pin-fin  channels  is  similar  for  the  1.0  mL/s  flow  rate  but  the  pressure 
drop  is  significantly  lower  for  the  taller  pin-fins  because  of  the  substantial  increase  in  cross- 
sectional  area  in  the  flow  direction.  The  surface  heat  transfer  coefficient  for  the  taller  pin-fins  is 
less  (attributed  to  the  decrease  in  flow  velocity)  and  these  base  heat  transfer  coefficient  results 
imply  that  the  increase  in  surface  area  counteracts  the  decrease  in  the  surface  heat  transfer 
coefficient  for  the  taller  pin-fins.  However,  the  substantial  decrease  in  the  pressure  drop  is 
especially  important  for  applications  where  pumping  power  must  be  minimized.  Overall,  these 
preliminary  results  indicate  that  taller  pin-fins  are  ideal  for  low  pumping  power  applications; 
however,  more  data  needs  to  be  gathered  before  drawing  any  definite  conclusions. 

As  shown  in  Figure  8,  while  the  experimental  heat  transfer  coefficient  results  are  generally  in  the 
range  of  the  correlation-derived  predictions,  the  average  percent  error,  at  greater  than  1 10%,  is 
unacceptably  high.  The  correlation  under  predicts  the  heat  transfer  coefficient  for  most  cases  and 
especially  at  higher  flow  qualities.  Also,  the  trends  contradict  one  another;  the  correlation 
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predicts  a  heat  transfer  coefficient  that  falls,  on  average,  with  increasing  flow  quality,  whereas 
the  experimental  results  demonstrate  a  heat  transfer  coefficient  that  steadily  increases  with  flow 
quality.  Accordingly,  some  additional  care  -  beyond  the  application  of  the  Reeser  et  al 
correlation  -  must  be  taken  in  designing  the  LD  micro  pin-fin  cooler  to  reflect  these 
experimental  trends. 


Figure  8:  Base  Heat  Transfer  Coefficient  vs.  Exit  Quality  for  Two-Phase  Micro  Pin-Fin 

Channel 

FC-72,  0.5  to  1.0  mL/s,  12  mm  x  12  mm  channel,  inline  pin-fin  configuration,  pin-fins  are  150 
pm  wide,  either  100  pm  or  500  pm  tall,  and  have  a  pitch  of  300  pm.  Dotted  line  data  is 
predicted  heat  transfer  coefficient  using  pin-fin  correlation  from  [2], 


Figure  9:  Pressure  Drop  vs.  Exit  Quality  for  Two-Phase  Micro  Pin-Fin  Channel 

FC-72,  0.5  to  1.0  mL/s,  12  mm  x  12  mm  channel,  inline  pin-fin  configuration,  pin-fins  are  150pm 
wide,  either  100pm  or  500pm  tall,  and  have  a  pitch  of 300pm. 
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2.2  Apparatus  of  Cryogenic  Microcoolers  for  High-Power  LD  Bar 
2.2.1  Final  Design  of  Micro-Pin  Fin  Manifold 

Armed  with  the  insights  obtained  from  the  correlation  study  and  FC-72  micro-pin  fin 
benchmarking  experiment,  a  final  liquid  nitrogen  cooler  design  was  produced.  A  cross-sectional 
and  exploded  view  of  the  final  design  iteration  is  shown  in  Figure  10  and  11,  respectively.  The 
cooler  consists  of  three  copper  parts:  an  upper  insert,  a  lower  insert,  and  a  central  manifold.  The 
micro  pin-fins  are  milled  out  of  the  upper  insert  and  positioned  under  the  LD  bar  or  test  heater. 
The  design  parameters  of  the  pin-fin  array  itself  are  identical  to  the  baseline  of  the  correlation 
study:  150  pm  wide,  300  pm  tall  pins  at  a  pitch  (streamwise  and  transverse)  of  300  pm. 
However,  the  footprint  of  the  array  has  been  reduced  to  10  mm  x  10  mm.  This  both  better 
matches  the  LD  bar  emitter  width  as  well  as  the  footprint  of  the  ceramic  heater  used  for  cooler 
characterization. 

The  central  piece  serves  as  a  plenum/manifold,  channeling  LN2  to  and  from  the  micro  pin-fin 
array.  Finally,  the  lower  insert  seals  off  the  lower  plenum  machined  in  the  underside  of  the 
central  manifold.  As  shown  in  the  exploded  view,  the  three  sections  of  the  cooler  are  bolted 
together  and  contact  surfaces  are  sealed  using  a  flexible  cryogenic  sealant. 

For  the  initial  evaluation  of  the  thermal  performance,  a  resistive  heater  is  placed  directly  above 
the  pin-fins  in  order  to  simulate  the  waste  heat  of  an  LD  bar.  The  modular  design  of  the  cooler 
facilitates  a  simple  exchange  of  different  pin-fin  designs  for  future  studies.  Following  the 
evaluation  of  the  cooler  thermal  performance,  the  LD  bar  is  then  bolted  to  the  non-wetted  side  of 
the  upper  insert,  using  a  thermal  interface  material  (TIM)  between  the  cooler  and  LD  bar  (as 
diagrammed  in  Figure  4). 


MOUNT  LD  BAR 
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HEATER 
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Figure  10:  Cross-Sectional  View  of  Cryogenic  Micro  Pin-Fin  Cooler 
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Figure  11:  Exploded  View  of  Cryogenic  Micro  Pin-Fin  Cooler 
2.2.2  Design  and  Final  Assembly  of  LN2  Flow  Loop  Apparatus 

The  cryogenic  flow  loop  apparatus  is  critical  in  both  regulating  the  supply  of  liquid  nitrogen  to 
the  pin-fin  array  and  in  collecting  experimental  data  during  both  the  test  heater  reference 
experiments  and  the  powered  LD  bar  tests.  While  the  flow  loop  underwent  several  iterations 
during  the  project,  the  basic  open-loop  design  philosophy  remained  constant. 

Shown  in  Figure  12  -  along  with  a  schematic  in  Figure  13  -  is  the  underlying  flow  loop  before 
installation  of  the  cooler  manifold  and  cryogel  insulation.  LN2  is  stored  and  supplied  via  a 
pressurized  tank,  where  the  tank’s  integrated  pressure  regulator  is  used  to  regulate  the  flow  of 
LN2  through  the  flow  loop  and  cooler.  The  LN2  flows  from  the  tank  into  a  four-way  junction, 
where  the  temperature  (Tinlet)  and  pressure  (Pinlet)  are  measured  with  an  Omega  E-Type 
thermocouple  and  cryogenic  pressure  transducer,  respectively.  The  LN2  then  proceeds  to  flow 
into  the  LD  cooler  (not  pictured,  installed  at  the  blue  end-caps  visible  in  Figure  12),  where  the 
flow  absorbs  heat  from  the  LD  bar  and  phase  changes  from  liquid  to  vapor.  The  two-phase 
mixture  then  exits  the  LD  cooler  and  flows  through  another  four- way  junction,  where  the 
temperature  (Toutlet)  and  pressure  (Poutlet)  are  measured  again.  The  pressure  and  temperature 
sensors  monitor  the  saturation  temperature,  bulk  fluid  temperature,  and  pressure  drop  across  the 
LD  two-phase  cooler. 

After  exiting  the  four-way  junction,  the  two-phase  flow  of  LN2  then  proceeds  through  a  liquid- 
to-liquid  heat  exchanger.  The  purpose  of  the  heat  exchanger  is  to  completely  evaporate  all  of  the 
excess  LN2  and  heat  the  nitrogen  gas  to  roughly  room  temperature,  so  that  the  mass  flow  rate  can 
be  measured  with  the  nitrogen  gas  mass  flow  meter.  The  LN2  is  heated  with  the  heat  exchanger 
using  a  constant  temperature  water  loop  maintained  at  25°C.  The  temperature  of  the  nitrogen  gas 

12 

Approved  for  public  release;  distribution  is  unlimited. 


is  measured  after  the  heat  exchanger  with  an  Omega  E-Type  thermocouple  to  ensure  complete 
evaporation  and  heating  of  the  LN2  to  standard  atmospheric  conditions. 

All  of  the  Swagelok  compression  fittings  and  tubing  are  made  of  316SS  and  the  apparatus  is 
insulated  using  Cryogel  Z  insulation.  Data  is  measured  and  logged  in  Lab  VIEW  using  a  National 
Instruments  NI  9214  (thermocouple)  and  NI  9205  (voltage)  cDAQ  module. 


Figure  12:  Assembled  LN2  Flow  Loop  Apparatus 
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2.2.3  LN2  Flow  Loop  Apparatus  Heat  Leakage  -  Cryogel  Z  Insulation 

Because  the  saturation  temperature  of  liquid  nitrogen  is  so  much  lower  than  room  temperature 
(-196°C  at  1  bar),  heat  leakage  from  the  surrounding  ambient  into  the  LN2  apparatus  must  be 
minimized.  Minimizing  the  heat  leakage  will,  in  turn,  reduce  the  flow  quality  at  the  inlet  by 
preventing  premature  (parasitic)  phase  change,  maximizing  the  amount  of  latent  heat  available 
for  the  LD  micro  pin-fin  cooler.  While  frequently  a  sealed  chamber  under  high  vacuum  can  be 
used  to  thermally  insulate  a  cryogenic  apparatus,  this  approach  is  very  expensive  and  complex.  A 
more  direct  approach  is  to  use  insulation  designed  for  cryogenic  applications,  such  as  Cryogel  Z 
(made  by  Aspen  Aerogels).  A  plot  of  thermal  conductivity  versus  temperature  for  Cryogel  Z 
insulation  is  shown  below  in  Figure  14,  where  it  can  be  observed  that  the  thermal  conductivity 
falls  with  decreasing  temperature  and  has  a  temperature-averaged  value  of  about  0.015  W/m-K, 
approximately  half  the  value  of  dry  air.  The  claimed  thermal  conductivity  value  was  also 
experimentally  confirmed  for  an  average  insulation  temperature  of  35  and  45°C.  The 
experimentally  determined  thermal  conductivity  was  only  slightly  lower  than  the  thermal 
conductivity  claimed  in  the  datasheet  and  within  the  limits  of  experimental  error. 

The  amount  of  heat  leakage  for  the  LN2  flow  loop  was  approximated  for  the  first  portion  of  the 
flow  loop,  consisting  of  all  the  components  leading  up  to  the  LD  cooler  inlet  and  the  LD 
bar/cooler  itself.  For  the  heat  leakage  approximation,  it’s  assumed  that  the  inner  wall  temperature 
(e.g.,  wall  temperature  of  stainless  steel  plumbing  and  fittings)  is  -196°C,  ambient  temperature  is 
25  °C,  and  the  dominate  mode  of  heat  transfer  is  natural  convection.  A  1-D  heat  transfer  analysis 
was  then  performed  using  the  known  geometry  of  the  plumbing,  fittings,  and  fixtures.  The  results 
from  the  heat  leakage  approximation  are  shown  in  Figure  15  with  a  plot  of  cooler  inlet  quality, 

X,  versus  flow  rate  of  LN2.  As  shown,  without  insulation  the  inlet  quality  is  very  high  at  lower 
flow  rates  (e.g.,  45%  at  a  flow  rate  of  1  mL/s),  indicating  that  the  heat  leakage  is  significant  and 
thus  reducing  the  amount  of  latent  heat  available  for  the  LD  cooler.  Conversely,  the  cooler  inlet 
quality  is  significantly  reduced  with  just  5  mm  of  Aerogel/Cryogel  Z  insulation  and  further 
reduced  with  10  mm  of  insulation.  With  10  mm  of  Cryogel  Z  insulation,  the  cooler  inlet  quality 
is  estimated  at  8%  for  a  flow  rate  of  1  mL/s  -  an  acceptable  coolant  expenditure  given  the 
remaining  latent  heat  available  in  the  remaining  liquid  delivered  to  the  pin-fin  array. 
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Temperature  (°F) 


Temperature  (°C) 


Figure  14:  Thermal  Conductivity  vs.  Temperature  for  Cryogel  Z  Insulation 

Source:  Aspen  Aerogels  datasheet. 
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Figure  15:  LD  Cooler  Inlet  Quality  versus  Flow  Rate  from  Heat  Leakage  Approximation 

for  LN2  Flow  Loop  Apparatus 


2.2.4  Preliminary  Testing  of  LN2  Flow  Loop  Apparatus 

Following  the  assembly  of  the  LN2  flow  loop,  a  cursory  two-phase  cooling  test  was  performed. 

In  order  to  quickly  evaluate  the  functionality  of  the  flow  loop  and  to  assess  whether  the  loop  is 
performing  as  designed,  a  Lytron  off-the-shelf  aluminum  cold  plate  was  attached  to  the  LN2  flow 
loop  in  place  of  the  pin-fin  cooler.  A  10  mm  x  10  mm  ceramic  resistive  heater  and  E-type 
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thermocouple  was  attached  to  the  cold  plate  in  order  to  evaluate  the  two-phase  thermal 
performance  of  the  cold  plate  configuration  with  LN2,  as  shown  in  Figure  16. 


Figure  16:  Lytron  Cold  Plate  Setup  used  to  Test  LN2  Flow  Loop  Apparatus 

The  initial  testing  of  the  flow  loop  was  satisfactory  in  that  it  didn’t  reveal  significant  problems 
with  the  functionality  of  the  flow  loop.  The  LN2  tank  provided  a  very  constant  source  of  pressure 
and  flow  -  with  very  little  fluctuations.  The  tank’s  valve  allowed  precise  control  over  the  flow 
rate.  Moreover,  the  pressure  transducers,  thermocouples,  and  mass  flow  meter  performed  as 
expected  and  the  water-nitrogen  heat  exchanger  effectively  evaporated  and  superheated  any 
remaining  LN2  to  room  temperature. 

A  plot  of  the  heat  transfer  coefficient  (HTC)  versus  calculated  change  in  flow  quality  (outlet 
quality  -  inlet  quality)  for  the  Lytron  cold  plate  with  LN2  is  shown  in  Figure  17.  This  HTC  is 
calculated  using  the  temperature  rise  of  the  plate-heater  interface  with  respect  to  the  inlet  fluid 
saturation  temperature,  with  the  footprint  of  the  ceramic  heater  used  as  the  normalization  area. 
No  insulation  was  used  for  this  initial  test  so  the  inlet  quality  is  not  zero  (or  near  zero),  hence  the 
x-  axis  is  given  as  a  change  in  flow  quality  and  not  an  absolute  value.  The  heat  transfer 
coefficient  profile  shown  in  Figure  17  is  for  a  flow  rate  of  3.0  mL/s  (2.42  g/s)  of  LN2,  an  inlet 
saturation  temperature  of  -184°C  (slightly  above  atmospheric  pressure),  pressure  drop  from  8.5 
to  14.9  kPa,  and  heat  flux  of  50  to  150  W/cm2.  As  shown,  the  HTC  rises  with  increasing  flow 
quality  and  heat  flux  up  to  a  change  in  flow  quality  of  around  10%,  where  it  peaks  at  a  value  of 
15,000  W/m2-K.  The  heat  transfer  coefficient  then  quickly  falls  off  for  changes  in  flow  quality 
above  10%,  reaching  a  minimum  of  6500  W/m2-K  at  a  change  of  flow  quality  of  14.2%.  This 
rapid  drop  in  the  heat  transfer  coefficient  is  likely  caused  by  dryout  of  the  heated  surface. 
Nonetheless,  the  profile  obtained  is  similar  in  shape  and  magnitude  to  the  thermal  data 
previously  published  by  Qi  et  al.  [3]  for  flow  boiling  of  LN2  in  mini  tubes. 
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This  same  basic  procedure  is  used  in  future  tests  to  evaluate  the  thermal  performance  of  the 
micro  pin-fin  cooler. 


Figure  17:  Heat  Transfer  Coefficient  versus  Change  in  Flow  Quality  (outlet  -  inlet)  for 

Lytron  Cold  Plate  with  LN2 

Q  =  3.0  mUs,  G  =  110  kg/m2-s,  q”  =  50-150  W/cm2,  Tsat.avg  =  -184°C,  A P  =  8.5-14.9  kPa,  and 

D  =  8mm  (copper  tube  with  U-bend). 

2.2.5  Installation  and  Initial  Testing  of  LN2  Pin-Fin  Cooler 

Following  the  assembly  of  the  flow  loop  and  preliminary  cold  plate  testing,  the  micro  pin-  fin 
cooler  was  installed  and  the  flow  loop  insulated.  The  installed  cooler  is  shown  in  Figure  1 8 
occupying  the  location  of  the  prior  Lytron  cold  plate.  As  with  the  cold  plate,  a  10  mm  x  10  mm 
ceramic  resistive  heater  is  attached  to  evaluate  the  two-phase  thermal  performance  of  the  pin-fin 
cooler.  Three  thermocouple  holes  provide  temperature  readings  at  the  inlet,  midpoint,  and  outlet 
of  the  pin-fin  array,  and  are  used  in  conjunction  with  dissipated  power  of  the  heater  to  determine 
the  pin- fin  cooler’s  performance. 
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Figure  18:  Assembled  LN2  Pin-Fin  Cooler 

For  the  two-phase  flow  process,  in  the  absence  of  an  ability  to  empirically  determine  the  input  or 
output  flow  quality,  the  change  in  quality  across  the  microcooler  is  the  only  determinable  quality 
parameter.  As  discussed,  in  order  to  provide  flow  to  the  pin-fin  cooler  with  near- zero  inlet 
quality,  aerogel  insulation  was  installed  on  the  flow  loop  and  manifold.  The  installed  cooler 
insulation  is  shown  in  Figure  19.  Further  aerogel  insulation  (not  shown  in  Figure  2)  was  added  to 
the  cryogenic  hose  supplying  LN2  from  the  tank. 
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Figure  19:  Flow  Loop  and  Manifold  Aerogel  Insulation 


2.2.6  Design  and  Purpose  of  Bypass/Blowdown  System 

A  final  flow  loop  modification  was  made  to  address  the  difficulty  of  bringing  the  liquid 
vaporization  front  into  and  past  the  pin-fin  array  with  the  available  22  psig  (237  kPa  abs) 
delivered  by  the  nitrogen  tank.  Because  the  system  startup  from  room  temperature  involves 
single-phase  gas  flow,  the  pressure  drop  required  for  a  given  startup  mass  flow  through  the  array 
are  much  higher  than  if  the  nitrogen  flow  was  in  liquid  or  two-phase  conditions.  This  is  in 
contrast  to  typical  two-phase  applications  where  the  saturation  temperature  is  above  room 
temperature,  where  startup  begins  with  zero-quality  liquid  and  transitions  to  two-phase  after  the 
thermal  load  to  the  cooler  is  applied.  Without  a  way  to  bypass  this  gas  flow  restriction,  the  initial 
cool-down  of  the  pin-fin  cooler  and  manifold  occurs  on  the  order  of  an  hour  or  more,  depending 
on  the  LN2  supply  tank  pressure. 


To  more  rapidly  bring  the  system  down  to  cryogenic  temperatures,  a  bypass  was  installed  prior 
to  the  pin-fin  manifold  inlet  temperature  and  pressure  fitting,  as  diagrammed  in  Figure  20.  This 
allows  a  much  increased  flow  of  LN2  (not  measured  but  easily  in  excess  of  15  mL/s),  rapidly 
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cooling  down  upstream  components.  In  addition  to  this,  the  exhaust  from  this  bypass  is  directed 
through  two  parallel  runs  of  14”  copper  tubing  clamped  on  either  side  of  the  copper  pin-fin 
manifold  before  venting  to  atmosphere  where  the  tubes  emerge  from  the  aerogel  insulation.  In 
this  way  the  pin-fin  array  itself  is  brought  down  to  cryogenic  temperatures  during  a  span  of 
several  minutes.  Upon  closing  the  bypass  valve,  the  liquid- vapor  front  is  now  past  the  pin-fin 
array.  Thus  the  fluid  passing  through  the  array  is  mostly  liquid,  dramatically  reducing  the 
pressure  needed  to  drive  a  desired  mass  flux.  With  the  heater  off  and  the  50  psig  tank,  the  flow 
obtained  through  the  loop  was  2.8  mL/s  of  liquid  -  or  2.25  g/s  -  with  a  pin-fin  manifold  inlet-to- 
outlet  pressure  drop  of  214  kPa  (with  a  further  55  kPa  being  dropped  from  outlet  to  ambient 
within  the  evaporative  heat  exchanger  and  outlet  tubing). 
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Figure  20:  Diagram  of  Bypass/Blowdown  Flow  Modification 

Opening  the  bypass  valve  directs  liquid  and  vapor  N2  through  tubing  mounted  on  outer  surface 

of  cooler  manifold. 


2.3  Performance  of  Cryogenic  Microcoolers 
2.3.4  Micro  Pin-Fin  Array  Initial  Two-Phase  Testing 

With  the  higher  pressure  tank  and  the  flow  loop  modification  allowing  adequate  liquid  delivery 
to  the  cooler,  a  set  of  initial  heater-based  tests  was  conducted.  From  the  outset,  an  encouraging 
observation  was  that  even  with  substantial  heater  power,  the  temperature  rise  measured  by  the 
embedded  thermocouples  remained  below  3  to  4°C,  indicating  the  feasibility  of  using  the  system 
to  maintain  an  LD  bar  at  cryogenic  temperatures.  Figure  21  quantifies  this  result  with  the  cooler 
manifold  conductance  for  three  stable  high-power  operating  points,  with  the  chip  heater  driving 
voltage  set  to  40,  60,  and  70  Volts.  The  chip  power  (voltage-current  product)  is  used  with  the 
measured  flow  rate  at  each  point  to  estimate  the  change  in  fluid  quality  just  as  with  the  cold  plate 
testing.  The  40,  60,  and  70  Volt  cases  correspond  to  the  7.5%,  19%,  and  30%  quality  points  in 
Figure  21  (chip  power  of  31,  70,  and  96  W,  respectively).  The  thermal  conductance  of  the  pin-fin 
array  is  defined  as  the  chip  power  divided  by  the  average  temperature  rise  of  the  three  embedded 
manifold  thermocouples  located  just  above  the  pin-fin  heated  wall  (see  Figure  23)  with  respect  to 
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the  zero  chip  power  case.  This  is  the  best  way  to  estimate  the  wall  temperature  rise  with  respect 
to  the  local  fluid  saturation  temperature  with  the  given  thermocouple  locations.  In  Section  3.4 
additional  thermocouple  locations  are  added  that  provide  a  more  accurate  instantaneous  wall 
superheat  measurement. 
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Figure  21:  Pin-Fin  Array  and  Manifold  Conductance  as  a  Function  of  Estimated  Fluid 

Quality  Change 

Inlet  pressure  is  fixed  at  370  kPa  absolute,  manifold  outlet  pressure  increases  with  quality  due  to 

increased  pressure  losses. 

Normalizing  this  array  conductance  by  the  1.14  cm2  footprint  of  the  array  (dimensions  shown  in 
Figure  22)  results  in  the  very  large  value  of  -250,000  W/m2-K  base  equivalent  heat  transfer 
coefficient.  However,  there  are  two  primary  uncertainties  involved  in  the  calculation  of  this 
value.  The  first  is  that  the  5  mm  upper  manifold  insert  that  lies  between  the  chip  heater  and  the 
pin-fin  array  operates  as  a  thermal  spreader,  allowing  a  portion  of  the  chip  power  to  conduct  to 
parts  of  the  copper  manifold  far  from  the  array  and  substantially  increasing  the  “wetted”  area  for 
heat  transfer  to  the  flowing  nitrogen.  While  this  is  obviously  of  benefit  in  keeping  the  heater  as 
cold  as  possible,  it  means  that  using  the  1.14  cm2  footprint  of  the  array  is  a  substantial 
underestimate  of  the  effective  cooling  area,  which  could,  perhaps,  be  as  much  as  twice  as  large. 
The  other  uncertainty  is  in  knowing  precisely  what  the  nitrogen  saturation  temperature  is  at  the 
array.  While  the  thermocouple  readings  in  the  unpowered  case  provide  a  good  estimate  of  the 
saturation  temperature,  powering  the  chip  alters  the  static  pressure  within  the  array,  changing  the 
saturation  temperature.  Allowing  for  a  ±1°C  uncertainty  in  this  temperature  contributes  a  ±30% 
uncertainty  in  the  heat  transfer  coefficient,  due  to  the  low  (3  to  4°C)  temperature  rise  observed. 

In  future  studies  addressing  these  uncertainties  can  be  done  on  two  fronts.  Uncertainty  in  Tsat  can 
be  reduced  by  embedding  thermocouples  below  the  pin-fin  array,  where  a  zero-flux  condition 
can  be  assumed  at  the  lower  wall.  To  accurately  predict  the  portion  of  chip  power  dissipated 
through  the  array  itself,  numerical  modeling  -  calibrated  with  the  available  temperatures  and 
used  to  determine  the  local  heat  fluxes  and  wall  temperatures  -  may  be  the  most  straightforward 


21 

Approved  for  public  release;  distribution  is  unlimited. 


approach.  Additional  thermocouple  wells  were  drilled  in  the  manifold  for  the  final  LD  bar 
testing.  These  thermocouple  locations  are  shown  in  Figure  23. 


Figure  22:  Central  Manifold  Block  Dimensions  (inches) 

The  footprint  of  the  array  and/or  micro  gap  are  computed  from  the  wetted  dimensions  of  the 
restriction  created  by  this  part,  0.4”  wide  x  0.443”  long,  with  0.013”  (300  pm)  gap  height. 
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Figure  23:  Thermocouple  Locations  on  the  Micro  Pin-Fin  Cooler  Manifold 

A  challenge  encountered  in  the  current  pin-fin  flow  loop  is  that  a  fixed  inlet  pressure  is  available 
from  the  LN2  tank,  rather  than  a  user-selectable  fixed  flow  rate  that  would  be  obtained  from  a 
positive  displacement  pump.  This  means  that  without  any  flow  control  within  the  loop  the  flow 
through  the  pin-fin  array  varies  as  chip  power  is  applied  due  to  an  increasing  pressure 
differential.  Plotted  in  Figure  24  is  the  nitrogen  mass  flowrate  as  a  function  of  applied  chip 
power. 

Heat  dissipated  from  the  chip  vaporizes  increasing  amounts  of  LN2,  leading  to  increased  pressure 
drop  from  both  the  expansion  and  acceleration  of  the  gas  phase.  At  80  volts  of  chip  power 
(126  W),  the  decrease  of  flow  is  sufficient  to  cause  a  runaway  increase  in  pressure  drop  as  falling 
mass  flow  leads  to  higher  exit  quality,  leading  to  additional  pressure  loss  which  eventually 
pushes  the  vaporization  front  upstream  of  the  pin-fin  array.  With  only  N2  vapor  as  the  coolant, 
the  heater  temperature  rises  continuously.  Fortunately,  unpowering  the  chip  heater  leads  to  a 
gradual  recovery  of  liquid  flow  in  a  matter  of  minutes,  indicating  that  there  is  sufficient  thermal 
mass  in  the  manifold  block  to  avoid  the  need  to  restart  the  bypass  and  blowdown  procedure.  An 
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automated  flow  control  method  would  be  able  to  prevent  this  flow  instability,  but  would  need  to 
down-  regulate  heater  power  if  a  prescribed  maximum  pressure  differential  is  reached. 
Fortunately  for  the  LD  bar  cooling  application,  126  W  is  well  over  the  expected  thermal 
dissipation  to  be  managed,  by  a  factor  of  2  or  more. 


Pin-Fin  Array  Flow  Vs  Chip  Power 
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Figure  24:  Fluid  Mass  Flow  Rate  as  a  Function  of  Applied  Chip  Heater  Power  (left  axis) 
and  Absolute  Pressure  at  Manifold  Inlet,  Outlet,  and  Ambient  (flow  loop  outlet  at  gas 

flowmeter)  (right  axis) 


2.3.2  Microgap  Two-Phase  Testing 

Two-phase  microgap  testing  was  undertaken  due  to  two  motivations.  First  it  offers  a  simple, 
low-pressure  alternative  to  the  pin-fin  array  that  could  be  tested  prior  to  the  solution  delivered  by 
the  bypass  loop  modification.  Second,  the  simplified  geometry  and  longer  history  of  use 
compared  to  micro  pin-fin  arrays  means  that  the  two-phase  correlations  that  will  be  compared  to 
experimental  results  are  simpler  and  more  mature  than  those  for  pin-fins.  Due  to  the  modular 
nature  of  the  pin-fin  manifold  design,  undertaking  microgap  experiments  required  only  the 
interchange  of  the  upper  and  lower  manifold  inserts.  The  pin-fin  array  is  “stored”  within  the 
lower  outlet  plenum,  protecting  it  from  damage  while  having  negligible  impact  on  the  nitrogen 
flow.  The  only  modification  made  was  the  addition  of  three  new  thermocouple  holes  to  capture 
the  heated  wall  temperature  of  the  new  microgap.  This  gap  has  the  same  height  as  the  pin-fin 
length  plus  clearance  height:  330  pm. 


Because  the  pressure  drop  across  the  microgap  is  substantially  less  than  that  for  the  pin-fin  array, 
no  startup  bypass  procedure  is  necessary.  Furthermore,  because  heater-induced  variation  in 
pressure  drop  is  also  decreased,  it  is  much  more  straightforward  to  regulate  mass  flowrate  using 
the  pressure  venting  valves  on  the  LN2  tank,  meaning  data  can  be  taken  at  constant  flow 
conditions  without  a  more  elaborate  flow  control  system.  As  a  preliminary  test  of  two-phase 
microgap  cooling,  results  were  collected  at  two  nominal  flow  rates  of  3.2  g/s  and  2.3  g/s  at 
different  chip  heater  powers  ranging  from  0  -  60  W.  The  results  are  presented  in  Figure  25 
through  Figure  28. 
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Cooler  Conductance  vs  Chip  Power 
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Figure  25:  Microgap  Cooler  Heat  Transfer  Coefficient  as  a  Function  of  Wall  Heat  Flux  for 

Two  Mass  Flowrates 
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Pressure  Drop  vs  Chip  Power 
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Figure  26:  Microgap  Cooler  Pressure  Drop  as  a  Function  of  Wall  Heat  Flux  for  Two  Mass 

Flowrates 
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Cooler  Conductance  vs  Quality 
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Figure  27:  Microgap  Cooler  Heat  Transfer  Coefficient  as  a  Function  of  Change  in  Vapor 

Quality  for  Two  Mass  Flowrates 
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Figure  28:  Microgap  Cooler  Pressure  Drop  as  a  Function  of  Change  in  Vapor  Quality  for 

Two  Mass  Flowrates 


In  these  tests,  the  cooler  conductance  is  defined  the  same  way  as  in  the  pin-fin  testing:  the 
average  of  the  temperature  rise  experienced  by  the  three  thermocouples  located  at  the  heated 
microgap  wall  with  respect  to  the  unpowered  condition.  If  the  measured  conductance  of  the 
microgap  cooler  were  to  be  normalized  by  the  same  1.14  cm2  as  the  pin-fin  array,  the  resulting 
equivalent  heat  transfer  coefficient  is  -100,000  W/m2-K.  As  with  the  pin-fin  array,  these 
equivalent  values  are  likely  higher  than  the  actual  local  values  since  thermal  spreading  within  the 
upper  insert  provides  a  parallel  path  to  parts  of  the  manifold  block  in  contact  with  (albeit  slower- 
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moving)  LN2  fluid.  This  spreading  path  is  likely  to  have  increased  participation  relative  to  the 
pin-  fin  application  due  to  the  lower  effective  heat  transfer  coefficient,  h,  of  the  microgap,  which 
leads  to  elevated  wall  and  heater  temperatures  compared  to  the  pin-fin  array.  Also  present  are  the 
same  sources  of  uncertainty  in  estimating  local  Tsat  within  the  gap,  although  the  lower  pressure 
variation  coupled  with  the  higher  temperature  differentials  lead  to  a  lower  preliminary 
uncertainty  of  ±20%.  The  same  measures  proposed  for  mitigating  pin-fin  Tsat  and  spreading 
uncertainty  can  be  applied  to  microgap  testing.  Despite  these  uncertainties,  one  indication  that 
appropriate  operation  is  being  obtained  is  that  the  ratio  of  pin-fin  to  microgap  cooler  conductance 
(and  by  extension  heat  transfer  coefficients)  is  between  2.5  -  2.8,  which  is  very  close  to  the  2.75 
time  area  enhancement  offered  by  the  pin-fin  array  over  the  microgap 

Regardless,  the  obtained  cooler  conductance  of  1 1  W/K  indicates  the  even  the  microgap  cooler  is 
certainly  capable  of  managing  a  substantial  amount  of  dissipated  power.  Identifying  the  thermal 
resistance  from  emitter-to-nitrogen:  Rld  is  0.8  K/W,  Ramch  estimated  at  0.24  K/W  for  silver 
grease,  Rcopper  is  400y.-,r_g.l-.w,iS„,iT,  =  0.09  K/W  and  Rcooier  is  0.09  K/W  results  in  a  total 
resistance  of  1.22  K/W,  of  which  the  microgap  cooler  is  less  than  10%.  This  is  a  particularly 
interesting  result:  until  the  other  resistances  between  the  emitter  and  nitrogen  are  addressed  (i.e. 
by  integrating  the  emitter  chip  and  submount  on  the  cooler  surface,  rather  than  relying  on  a 
separate  conduction  bar),  a  microgap  may  be  sufficient  is  preventing  the  convection  resistance 
from  becoming  the  thermal  bottleneck. 

One  drawback,  however,  is  that  although  the  microgap  is  simpler  and  presents  a  much  lower 
pressure  drop  than  the  pin  fins,  its  performance  is  measured  at  generally  higher  mass  flow  rates 
than  those  required  for  the  pin-fin  array.  Higher  flow  means  the  exit  quality  is  lower  than  the  in 
the  case  of  the  pin-fin  array,  indicating  that  a  smaller  portion  of  the  available  coolant  is  being 
used  for  thermal  management.  This  is  a  particular  issue  considering  the  open  loop  configuration 
of  the  LN2  flow  apparatus,  and  represents  a  large  loss  in  overall  efficiency  unless  the  un¬ 
vaporized  coolant  can  be  put  to  further  use  downstream  of  the  microgap  in  any  eventual 
application.  Redesigning  the  flow  apparatus  to  deliver  lower  flow  rate  and  thus  higher  exit 
quality  will  provide  insight  on  whether  the  microgap  maintains  similar  performance  compared  to 
high-flow  conditions. 

One  final  consideration  identified  in  the  preliminary  testing  can  be  seen  in  the  pressure  curves 
plotted  in  Figure  26  and  Figure  28.  The  tests  were  conducted  with  the  high  flow  first,  moving 
from  high  to  low  chip  power.  This  was  done  because  flow  control  by  tank  venting  is  relatively 
responsive  when  lowering  tank  (and  thus  inlet)  pressure,  but  requires  time  for  pressure  to  build 
back  up  to  regain  high  pressures.  For  the  3.2  g/s  flow,  as  chip  power  was  reduced  the  vapor 
quality  exiting  the  cooler  decreased,  meaning  more  LN2  need  to  be  evaporated  within  the  water- 
stabilized  heat  exchanger  before  exhausting  to  the  gas  flowmeter.  For  the  lowest  power  tests,  this 
eventually  overwhelmed  the  heater  responsible  for  water  temperature  control:  the  reservoir  water 
temperature  was  less  than  10°C,  and  icing  was  visible  on  the  outer  surfaces  of  the  heat 
exchanger.  N2  exhaust  temperature  prior  to  the  flow  meter  decreased  to  150  K  by  the  10  W 
operating  point,  well  outside  the  accurate  range  of  the  gas  flow  meter.  It  is  likely  mass  flow  was 
substantially  higher  than  measured  3.2  g/s  (3.9  mL/s),  as  evidenced  by  the  increasing  pressure  as 
chip  power  was  reduced.  The  flow  loop  was  shut  down  and  water  temperature  allowed  to  recover 
before  proceeding  with  the  2.3  g/s  testing,  which  displayed  relatively  constant  pressure  drop. 
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Exhaust  temperature  for  the  2.3  g/s  testing  remained  above  280  K  for  the  entire  test  run.  The 
lesson  is  that  for  high  LN2  flow  testing,  additional  care  is  necessary  in  monitoring  the  auxiliary 
components  of  the  flow  loop. 

2.4  Application  of  Cryogenic  Microcoolers  for  High-Power  LD  Bar 
2.4.2  LN2  Cooled  Laser  Diode  Bar  Experimental  Setup 

The  laser  diode  bar  was  fixed  to  the  micro-pin  fin  cooler  following  the  initial  ceramic  heater 
cooler  characterization  tests.  Shown  in  Figure  29  is  the  DILAS  conduction  cooled  laser  diode 
bar,  along  with  the  locations  where  additional  thermocouples  were  placed  for  measurements 
during  the  experiments. 
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Figure  29:  DILAS  Conduction  Cooled  Laser  Diode  Bar 

The  characteristics  of  this  LD  bar  are  as  follows:  19  identical  GaAs  emitters,  each  with  a  width 
of  100  pm  with  a  center-to-center  spacing  of  500  pm  between  emitters,  yielding  a  fill-factor  of 
20%.  The  length  of  the  chip  is  2000  pm,  with  an  active  epitaxial  layer  thickness  of  4.6  pm  grown 
on  a  chip  substrate  thickness  of  107  pm.  This  chip  is  attached  to  positive  copper  electrode  with 
indium  solder  with  a  layer  thickness  of  10  pm.  The  maximum  optical  power  is  65  W  at  62.5  A, 
with  a  center  wavelength  of  976  nm  +/-  10  nm  and  an  overall  efficiency  of  63.4%.  The  spec 
sheet  thermal  resistance  from  the  junction  to  the  copper  block  is  0.8  K/W. 

Figure  30  shows  the  LD  bar  mounted  on  top  of  the  micro-pin  fin  cooler  using  the  four  mounting 
bolts,  situating  the  emitter  in  the  vicinity  of  the  micro-pin  fin  array,  and  directing  the  beam 
towards  the  optical  power  sensor  (not  shown).  The  ends  of  copper  bypass  exhaust  tubes  were 
then  angled  outward  to  avoid  exhausting  directly  onto  the  power  meter. 
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Figure  30:  Nitrogen  Flow  Loop  with  DILAS  LD  Bar  Mounted  on  Micro-Pin  Fin  Cooler 

A  single  test  cycle  using  the  LD  bar  consisted  of  a  room  temperature  startup  of  the  LD  bar 
(powered  at  60  A),  shortly  followed  by  nitrogen  flow  through  both  the  pin  fin  array  and  the 
bypass.  This  was  done  to  reach  a  near-emitter  temperature  of  46°C  before  cooling  progressively 
to  cryogenic  temperatures.  46°C  is  the  measured  near-emitter  temperature  during  a  reference  test 
using  a  thermoelectric  temperature-controlled  stage  held  at  25  °C  and  a  LD  current  of  60  A, 
which  is  the  standard  operating  conditions  specified  in  the  LD  bar  data  sheet.  The  flow  of 
nitrogen  through  both  the  bypass  and  pin  fin  array  is  adjusted  to  slowly  bring  the  cooler  and  LD 
bar  temperature  down  to  cryogenic  temperatures  and  allow  two-phase  flow  to  be  established 
fully  in  the  pin  fin  array.  At  this  point  the  temperature  at  the  interface  between  the  cooler  and  the 
LD  bar  is  about  -180°C.  The  cooldown  -  which  occurs  over  10  to  15  minutes  -  is  assumed  to  be 
slow  enough  that  the  effects  of  transient  temperature  change  on  the  diode  emitter  can  be 
neglected.  The  LD  bar  power  and  efficiency  measurements  are  thus  treated  as  quasi-steady  state. 

Figure  31  illustrates  some  of  the  considerations  encountered  during  the  cryogenic  tests.  In  order 
to  provide  access  for  the  LD  bar  power  input  and  allow  a  path  for  the  beam  output,  the  aerogel 
insulation  had  to  be  rolled  back  from  the  end  of  the  cooler.  As  the  cooler  reaches  subzero 
temperatures  ice  begins  to  build  up  on  exposed  surfaces,  including  the  LD  bar  and  its  emitting 
surface.  Preliminary  cryogenic  tests  revealed  that  once  this  emitter  ice  buildup  occurs  measured 
optical  power  drops  markedly.  To  combat  this  a  nozzle  directing  nitrogen  gas  towards  the  emitter 
in  a  shield  plume  was  installed  above  the  cooler.  This  nitrogen  was  sourced  from  the  head  space 
in  the  LN2  supply  tank,  and  so  was  also  at  cryogenic  temperatures.  In  this  way  the  shield  gas 
attempts  to  reduce  the  impact  of  the  ambient  temperature  and  humidity  on  cooler  operation. 
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2.4.3  LN2  Cooled  Laser  Diode  Bar  Experiment  Results 


Figure  32  shows  experimental  results  for  two  separate  cryogenic  cycles  on  an  LD  bar  run  at  60  A 
emitter  current.  Measured  power  values  are  plotted  as  a  function  of  the  cooler/conduction  block 
interface  temperature,  providing  a  comparison  to  measurements  made  using  the  25  °C  thermal 
stage.  Cycles  proceeded  from  room  temperature  to  cryogenic  temperature,  i.e.  from  right  to  left 
on  the  plots  in  Figure  32. 


As  temperature  decreased,  the  input  electrical  power  rose  as  the  diode  forward  bias  voltage 
increased  (from  1.55V  at  25°C  to  1.83V  at  -180°C).  For  the  first  portion  of  each  cycle,  the 
measured  output  optical  power  increased  with  decreasing  temperature  up  until  an  interface 
temperature  between  -75  to  -100°C,  where  it  peaks  at  68.9W.  Measured  output  power  then  fell  as 
temperatures  dropped  further,  even  falling  below  levels  measured  at  room  temperature.  These 
reductions  in  measured  power  coincided  with  progressive  ice  film  accumulation  that  occurred 
despite  the  nitrogen  shield  gas.  The  interpretation  is  that  ice  buildup  on  the  emitter  surface 
occludes  and/or  scatters  the  emitted  beam,  and  moreover  that  differences  in  the  icing  process 
accounts  for  the  discrepancy  between  data  from  individual  tests.  This  beam  interference 
manifests  as  a  reduction  in  the  effective  power  efficiency  of  the  LD  bar  in  terms  of  optical  power 
delivered  to  the  power  meter,  as  seen  in  the  second  plot  of  Figure  32.  The  peak  effective  power 
efficiency  in  each  run  was  68.9%  and  68.4%,  respectively,  occurring  at  an  interface  temperature 
of -50°C. 


Figure  31:  Cooler  and  LD  Bar  at  the  End  of  a  Test  Cycle 
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Figure  32:  Test  Data  for  Two  Cryogenic  Cycles  on  an  unused  LD  Bar 


The  discrepancy  between  the  electrical  input  power  and  the  measured  optical  power  is  recorded 
as  the  residual  power  not  transmitted  to  the  optical  target.  At  temperatures  above  -50°C  this 
power  is  identified  as  the  LD  bar  thermal  dissipation,  but  after  an  ice  film  is  generated  it  remains 
an  open  question  if  the  entirety  of  the  residual  power  goes  to  heating  the  LD  bar  or  if  some  is 
scattered  and  unaccounted  optical  power. 


Accounting  for  this  power  in  the  thermal  operating  characteristics  of  the  cooler  is  complicated  by 
the  unknown  heat  transfer  facilitated  by  the  shield  gas  jet,  the  bypass  tubing,  and  ambient.  An 
example  of  this  is  cooler  operation  at  -25°C  where  an  energy  budget  is  presented  in  Figure  33, 
before  ice  film  generation  and  when  the  flow  in  the  pin  fin  array  is  entirely  single-phase  gas. 
Because  the  nitrogen  at  the  cooler  inlet  and  outlet  is  superheated  gas,  calculating  the  net  heat 
removal  is  a  matter  of  computing  the  net  enthalpy  flow,  while  the  absence  of  ice  means  the  30W 
LD  bar  residual  power  can  be  completely  ascribed  to  thermal  dissipation.  Comparing  these 
energy  rates  to  the  rate  at  which  the  cooler/LD  bar  mass  cools  down  highlights  a  40  -  50  W 
discrepancy,  which  is  attributed  to  the  cooling  effect  of  the  cold  shield  gas  jet  and  the  nitrogen 
flowing  through  the  bypass.  As  temperatures  continue  to  fall,  this  cooling  effect  will  decrease  as 
the  exposed  surfaces  of  the  cooler/LD  bar  become  closer  to  the  gas  jet  and  bypass  temperatures 
and  further  from  the  ambient  temperature. 
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Figure  33:  Energy  Budget  for  LD  Bar  Cooler  Operation  at  -25°C  Interface  Temperature 
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Below  -25  °C  interface  temperature  this  calculation  is  made  more  complicated  as  the  nitrogen  at 
the  cooler  inlet  reaches  saturation.  As  the  inlet  flow  quality  is  unknown,  the  internal  heat 
removal  rate  cannot  be  estimated.  However,  once  the  cooler  outlet  also  reaches  saturation  (at  an 
interface  temperature  of  -180°C),  the  flow  in  cooler  is  fully  two-phase,  picks  up  to  a  mass  flow 
rate  of  -2  g/s,  and  reaches  steady  state.  At  this  point  the  cooler  operation  can  be  compared  to  the 
tests  performed  with  the  ceramic  chip  heater  where  the  internal  cooler  conductance  was 
measured  to  be  about  30  W/K.  The  measured  wall  superheat  at  this  point  is  1.5+0.5K, 
suggesting  that  the  cooler  is  removing  between  30  and  60  W.  The  residual  LD  power  is 
between  50  and  60  W,  indicating  that  on  the  order  of  10W  is  unaccounted  for.  This  power  may 
have  been  removed  by  the  bypass  and  shield  jet,  or  as  scattered  optical  power  not  captured  by  the 
optical  power  meter. 
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3.  PART  II:  POF-BASED  RELIABILITY  ASSESSMENT  OF 
LD/MICROCOOLER  SUBASSEMBLY  SUBJECTED  TO  CRYOGENIC 

OPERATING  CONDITIONS 

3.1  PoF  Model  for  LD/Microcooler  Subassembly 

The  concept  of  a  hierarchical  PoF-based  reliability  assessment  model  is  proposed.  The  model  is 
shown  in  Figure  34  for  the  LD  bar  cooled  by  the  cryogenic  pin-fin  cooler.  The  LD  bar  is 
eventually  used  as  an  optical  power  source  for  a  fiber  laser.  The  failure  of  the  fiber  laser  can  be 
separated  by  catastrophic  optical  mirror  damage  from  the  LD  bar  and  gradual  degradation  of  the 
optical  power.  For  the  gradual  degradation,  there  is  no  standard  failure  criterion.  We  have 
decided  the  failure  criterion  for  the  gradual  degradation  when  the  optical  power  is  less  than  80% 
of  the  initial  optical  power. 

The  gradual  degradation  of  the  optical  power  governed  by  the  pumping  efficiency.  The  pumping 
can  be  affected  by  optical  reliability  and  LD  bar  reliability.  The  optical  reliability  is  related  to 
degradation  of  fiber  (lasing  medium)  and  coupling  efficiency.  The  LD  bar  reliability  is  governed 
by  optical  power  change  from  chip  degradation,  center  wavelength  shift,  spectrum  width  change, 
and  beam  divergence  and  near-field  linearity.  The  beam  divergence  and  near-field  linearity  also 
affect  the  coupling  efficiency  in  the  optical  reliability. 

The  LD  bar  reliability  is  mainly  affected  by  the  junction  temperature  distribution,  which  is 
governed  by  the  thermal  reliability.  When  the  forward  current,  inlet  coolant  temperature,  and 
flow  rate  are  determined  from  the  operation  input,  the  LD  bar  have  the  initial  junction 
temperature  distribution.  The  thermal  reliability  affects  the  junction  temperature  distribution.  The 
thermal  reliability  can  be  separated  by  thermal  conduction  resistance  of  the  die  attach  and  the 
heat  transfer  coefficient  of  the  cooler.  The  thermal  conduction  resistance  of  the  die  attach  can  be 
increased  by  voids  from  the  electromigration  or  cracks  from  the  thermal  fatigue.  The  heat 
transfer  coefficient  of  the  cooler  can  be  degraded  by  erosion,  corrosion,  clogging  in  the  pin-fin 
channel. 
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Figure  34:  Hierarchical  Life  Prediction  Model  for  the  Cryogenic  Pin-Fin  Cooled  LD  Bar 


In  this  study,  the  crack  propagation  in  the  die  attach  due  to  the  thermal  fatigue  and  its  effect  on 
the  SPD)  change  are  demonstrated  for  the  PoF-based  reliability  assessment.  First,  the  junction 
temperature  distribution  measurement  is  essential  as  the  initial  step.  A  hybrid 
experimental/numerical  method  is  proposed  to  predict  the  junction  temperature  distribution  of  a 
high  power  LD  bar  (Section  3.2).  In  Section  3.3,  a  novel  method  is  proposed  to  predict  the  SPDs 
of  individual  emitters  by  deconvoluting  the  SPD  of  an  LD  array.  Eventually,  this  method  can  be 
utilized  to  predict  the  SPD  change  from  the  crack  propagation  of  the  LD  bar.  In  Section  3.4,  the 
crack  propagation  in  the  die  attach  due  to  the  thermal  fatigue  and  its  effect  on  the  spectral  power 
distribution  change  are  demonstrated  for  the  PoF-based  reliability  assessment. 

3.2  Method  for  Predicting  Junction  Temperature  Distribution  in  High-Power  LD  Bar 

A  hybrid  experimental/numerical  method  is  proposed  for  predicting  the  junction  temperature 
distribution  in  a  high  power  LD  bar  with  multiple  emitters.  A  commercial  water-cooled  LD  bar 
with  multiple  emitters  is  used  to  illustrate  and  validate  the  proposed  method.  A  unique 
experimental  setup  is  developed  and  implemented  first  to  measure  the  average  junction 
temperatures  of  the  LD  bar  emitters.  After  measuring  the  heat  dissipation  of  the  LD  bar,  the 
effective  heat  transfer  coefficient  of  the  cooling  system  is  determined  inversely  from  the 
numerical  simulation  using  the  measured  average  junction  temperature  and  the  heat  dissipation. 
The  characterized  heat  dissipation  and  effective  heat  transfer  coefficient  are  used  to  predict  the 
junction  temperature  distribution  over  the  LD  bar  numerically  under  high  operating  currents.  The 
results  are  presented  in  conjunction  with  the  wall-plug  efficiency  and  the  center  wavelength  shift. 
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3.2.1  Introduction 


As  higher  optical  power  is  demanded  for  advanced  applications,  more  closely-spaced  emitters 
with  higher  forward  current  are  used  in  LD  bars.  As  a  result,  the  junction  temperature  from  the 
center  to  the  edge  emitters  may  have  large  variations,  which  makes  the  center  wavelength  and 
wall-plug  efficiency  of  each  emitter  different  from  each  other. 

Several  junction  temperature  measurement  methods  for  low  power  LDs  or  light  emitting  diodes 
(LEDs)  have  been  proposed,  including  techniques  based  on  measurement  of  the  thermal 
resistance  [4],  wavelength-shift  [5,  6],  optical  power  output  [5],  and  forward-voltage  [5,  7-14], 
These  methods  are  applicable  only  when  the  junction  temperature  is  uniform.  Micro-Raman 
spectroscopy  [15-17]  can  be  used  to  measure  the  junction  temperature  distribution  by  measuring 
multiple  local  temperatures.  In  practice,  it  requires  a  complicated  experimental  setup  and  has 
limited  accuracy  (10  to  20°C)  [15-17], 

3.2.2  Laser  Diode  System 

This  section  is  devoted  to  the  description  and  the  electrical  characteristics  of  a  commercial  water- 
cooled  LD  bar  tested  in  the  study. 

3.2.2. 1  LD  Bar  Description 

The  commercial  LD  bar  system  (E11.4N-940.10-150C-S013.1:  DILAS)  is  shown  in  Figure  35a. 
The  LD  bar  consists  of  23  identical  GaAs  emitters.  The  fill-factor  is  50%;  each  emitter  is  200  pm 
wide  and  has  a  pitch  of  400  pm.  The  maximum  optical  power  at  160  A  is  160  W  with  a  center 
wavelength  of  930  nm. 

The  close-up  view  of  the  side  of  the  LD  bar  is  shown  in  Figure  35b.  The  GaAs  chip  (the 
epi-down  configuration)  is  mounted  on  a  CuW  submount  using  AuSn  die  attach.  The  submount 
made  of  CuW  (coefficient  of  thermal  expansion  (CTE):  6.5  ppm/°C)  is  placed  between  the  GaAs 
chip  (CTE:  6.4  ppm/°C)  and  the  water-cooled  microchannel  made  of  Cu  (CTE:  16.6  ppm/°C)  for 
the  stress-relieving  buffer  layer  to  reduce  the  thermal  stress  attributed  to  the  mismatch  in  the 
CTE  between  them  as  well  as  for  the  heat  spreader  [18-20].  The  specified  thermal  resistance 
from  the  junction  to  water  inlet  temperature  is  approximately  0.3  K/W  [21],  The  internal 
structures  of  the  commercial  microchannel  and  the  interfacial  resistance  are  not  available,  and 
thus  the  effective  water  heat  transfer  coefficient  for  this  commercial  microcooler  cannot  be 
determined. 
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Figure  35:  (a)  LD  Bar  with  Water-Cooled  MicroChannel  [21]  and  (b)  Side  View  of  the 

LD  Bar 


3.2.2.2  Calibration  Curve 

It  has  been  known  that  a  negative  linear  relationship  exists  between  the  junction  temperature  and 
the  forward  voltage  of  a  laser  diode  [7-14],  The  junction  temperature  at  the  operating  current 
can  be  measured  using  this  relationship  (known  as  the  “calibration  curve”),  and  it  is  called  the 
forward  voltage  method  [7-14], 
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The  “calibration  curve”  is  obtained  using  a  probe  current  much  lower  than  the  operating  current. 
If  the  probe  current  is  too  low,  the  forward  voltage  loses  the  negative  linear  relationship  at  high 
junction  temperatures  due  to  the  leakage  current  effect  [22].  On  the  other  hand,  if  it  is  too  high, 
the  loss  of  linearity  occurs  at  low  junction  temperatures  due  to  the  internal  series  resistance  [23], 
In  addition,  the  probe  current  should  be  as  low  as  possible  to  avoid  any  undesired  junction 
temperature  increase  while  obtaining  the  calibration  curve.  Every  LD  has  somewhat  different 
electrical  characteristics.  Thus,  it  is  important  to  determine  the  lowest  probe  current  that  provides 
the  desired  linearity  [14]. 

The  LD  bar  was  placed  inside  a  convection  oven  (EC1A:  Sun  Electronics  Systems)  and  the 
forward  voltage  was  measured  at  25,  35,  45,  55,  and  65°C  (with  an  accuracy  of  ±0.1  °C)  by  a  data 
acquisition  module  (DAQ:  USB-6212:  National  Instruments)  with  a  16-bit  resolution.  The 
maximum  operating  junction  temperature  was  estimated,  based  on  the  thermal  resistance  of 
0.3  K/W.  the  inlet  water  temperature  of  20°C,  and  the  measured  maximum  heat  dissipation  of 
84.7  W  (this  will  be  explained  further  in  Section  3.2.5. 1),  to  be  45.4°C.  The  calibration  curve 
measurement  was  repeated  at  various  probe  current  values  (from  20  mA  to  160  mA  with  an 
interval  of  20  mA). 

The  results  are  shown  in  Figure  36,  displaying  the  expected  linear  relationship  between  forward 
voltage  and  junction  temperature.  The  small  deviations  from  linearity  in  voltage  and  temperature 
at  65°C  are  summarized  in  Table  1.  The  junction  temperature  error  gradually  decreased  as  the 
probe  current  increased,  and  remained  virtually  the  same  after  120  mA.  The  probe  current  of 
120  mA  generated  a  heat  dissipation  of  only  142  mW  at  25°C  (6.2  mW  per  emitter),  which  was 
negligible  compared  to  the  heat  dissipation  produced  by  the  operating  current.  Thus,  the 
calibration  curve  obtained  at  120  mA  was  selected  for  junction  temperature  measurement.  The 
slope  and  the  y-intercept  of  the  calibration  curve  were  -1.21  mV/K  and  1.2104  V,  respectively. 


Figure  36:  Forward  Voltage  as  a  Function  of  Junction  Temperature 
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Table  1.  Junction  Temperature  Error  at  65°C  under  different  Probe  Currents 


If  [mA] 

Deviation  of  Vf  (mV) 

Error  in  Tj  [°C] 

160 

0.17 

0.1 

140 

0.16 

0.1 

120 

0.18 

0.1 

100 

0.41 

0.3 

80 

0.78 

0.6 

60 

0.80 

0.7 

40 

0.80 

0.7 

20 

1.60 

1.3 

3.2.2.3  Electrical  Resistance  of  Single  Emitter 


The  calibration  curve  is  obtained  when  all  the  emitters  have  the  same  temperature.  The  emitters 
of  the  LD  bar  are  connected  in  parallel,  and  thus  the  electrical  resistance  of  a  single  emitter 
(assuming  that  all  emitters  are  identical)  can  be  determined  simply  by 


R(T)  =  Nx 


(1) 


where  R  is  the  electrical  resistances  of  the  single  emitter  [Q];  Vf  is  the  forward  voltage  of  the  LD 
bar  [V]  under  the  probe  current,  Iprobe  [A];  and  N  is  the  number  of  the  emitter  ( N  =  23). 


The  results  obtained  for  the  LD  bar  at  the  probe  current  of  Iprobe  =120  mA  are  shown  in 
Figure  37.  The  electrical  resistance  decreased  with  the  temperature;  the  change  in  resistance  was 
only  4%  (from  226.2  Q  to  216.9  Q)  over  the  temperature  range  from  25°C  to  65°C. 
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Figure  37:  Electrical  Resistance  of  the  Single  Emitter  as  a  Function  of  Junction 

Temperature  at  Iprobe  =  120  mA 
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3.2.3  Junction  Temperature  Measurement 

The  junction  temperature  at  the  operating  current  can  be  measured  by  switching  the  operating 
current  to  the  probe  current  [7-14].  As  discussed  in  Ref.  [24],  the  forward  voltage  shows  the 
combined  behavior  of  RC  delay  and  thermal  delay  during  the  switching  time.  The  RC  delay  is 
attributed  to  the  resistance  of  a  LD  and  the  capacitance  of  a  current  source. 

The  LD  bar  is  operated  at  very  high  forward  currents.  A  power  supply  that  drives  high  currents 
typically  has  large  capacitance,  which  can  cause  the  large  RC  delay,  and  the  transient  junction 
temperature  behavior  of  the  CuW  submount  cannot  be  documented.  Fast  switching  circuits  with 
two  separate  power  supplies  have  been  utilized  to  reduce  the  delay  [10,  11].  This  scheme  is 
adopted  for  the  current  study. 

3.2.3.1  Test  Setup 

A  test  apparatus  to  minimize  the  RC  delay  is  illustrated  schematically  in  Figure  38.  The 
operating  current  source  (LDX-36125-12:  ILX  Lightwave)  applies  the  operating  current  with  a 
nominal  accuracy  of  ±  (0.1%  +  120  mA).  The  probe  current  source  (2401:  Keithley  Instruments) 
applies  the  probe  current  with  a  nominal  accuracy  of  ±  (0.066%  +  20  pA).  The  two  power 
supplies  are  connected  in  parallel. 
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Figure  38:  Schematic  Illustration  of  Junction  Temperature  Measurement  Setup 


An  N-channel  metal-oxide-semiconductor  field-effect  transistor  (MOSFET)  (IRL7833PBF: 
International  Rectifier  Corporation)  serves  as  a  switch  for  the  operating  current  source  [10]. 
The  probe  current  flows  from  the  source  (S)  to  the  drain  (D)  of  the  N-channel  MOSFET,  even 
when  the  switch  is  off  to  block  the  operating  current.  A  diode  (150EBU02-ND:  Vishay 
Intertechnology)  is  inserted  between  the  operating  current  source  and  the  MOSFET  to  prevent 
this  undesired  flow.  The  MOSFET  and  the  diode  are  mounted  on  a  heat  sink  to  dissipate  heat  at 
high  operating  currents. 
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The  chiller  (ISOTEMP  1 1 15V/60HZ  PD-1:  Fisher  Scientific)  regulates  the  water  inlet 
temperature  with  a  temperature  stability  of  ±0.1  °C,  and  the  flow  regulator  (FLDW321 1G: 
OMEGA  Engineering)  controls  the  flow  rate.  The  optical  power  sensor  (USB-PM- 150-50: 
Coherent  Laser  Group)  measures  the  optical  power  with  a  nominal  accuracy  of  ±2.7%  in  the 
operating  range.  The  current  sources,  the  DAQs,  and  the  optical  power  sensor  are  integrated  into 
a  Lab  VIEW  program. 

In  the  actual  measurements,  the  chiller  is  set  to  produce  the  inlet  water  temperature  of  20°C  and 
the  flow  rate  of  16  L/h.  The  pressure  drop  from  the  water-cooled  microchannel  cooler  is  34  psi. 
The  forward  voltage  measurement  is  conducted  only  with  the  probe  current  applied  to  the  LD 
bar.  As  an  example,  to  apply  an  operating  current  of  80  A  to  the  LD  bar,  the  operating  current 
source  and  the  probe  current  source  applies  79.88  A  and  120  mA,  respectively,  with  the 
MOSFET  switch  “on”.  When  the  optical  power  and  the  forward  voltage  of  the  LD  reaches  the 
steady-state  condition,  the  switch  is  turned  off  to  block  the  flow  of  the  operating  current.  The 
data  acquisition  module  2  (DAQ2:  USB-6212:  National  Instruments)  supplies  the  gate  voltage  to 
the  MOSFET  switch  and  the  data  acquisition  module  1  (DAQ1:  USB-6212:  National 
Instruments)  measures  the  forward  voltage  of  the  LD  with  16-bits  resolution  with  the  maximum 
sampling  rate  of  400  kS/s  continuously  during  the  transient  period. 

3.2.3.2  Average  Junction  Temperature 

As  mentioned  earlier,  the  junction  temperatures  of  emitters  can  have  large  variations  at  high 
operating  currents.  However,  only  a  single  value  for  the  entire  LD  bar  can  be  obtained  from  this 
setup.  The  following  investigation  is  conducted  to  define  the  physical  meaning  of  the  measured 
value. 

The  emitters  are  connected  in  parallel,  and  thus  the  electrical  resistance  of  the  LD  bar  can  be 
expressed  as: 


1 


(2) 


„  <7)  = 


where  Rbar  is  the  electrical  resistance  of  the  bar  [Q];  Ri  and  Ti  are  the  electrical  resistance  [Q] 
and  the  temperature  of  the  emitter  [°C],  respectively. 

Let’s  consider  a  case  where  the  junction  temperature  of  the  LD  bar  increases  linearly  from  the 
edge  to  the  center  (AT  =  Tcenter  -  Tedge).  This  simple  linear  variation  is  analyzed  to  illustrate  the 
physical  meaning  of  the  measured  value.  The  average  junction  temperature  of  this  case  can  be 
expressed  as 
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Then,  the  true  forward  voltage  and  the  forward  voltage  estimated  based  on  the  average 
temperature  can  be  expressed  as: 


Vr(T)  =  ^(T)I^ 


V™(T)=±R(T™)Iprai 


(3) 


J 7*™  V™  and  —R(Tpr) 

f  is  the  true  forward  voltage  of  the  LD  bar  [V];  N  '  are  the  forward  voltage 

[V]  and  the  electrical  resistances  [Q]  of  the  LD  bar  at  the  average  junction  temperature, 
respectively.  It  is  to  be  noted  that  the  true  forward  voltage  and  the  forward  voltage  based  on  the 
average  temperature  are  different  because  the  emitters  are  connected  in  parallel. 


The  difference  between  these  two  values  provides  an  estimate  about  how  the  true  junction 
temperature  deviates  from  the  average  junction  temperature  of  the  bar  by  dividing  slope  of  the 
calibration  curve  [V/K].  The  deviation  can  be  defined  as: 


ST. 


T/mtnjs  Trr-avtf’ 

jf  ~  V 

1.21  L0“3 


(4) 


where  STj  is  deviation  of  the  true  junction  temperature  from  the  average  junction  temperature  of 
LD  bar  [°C];  1.21-10'3  is  the  slope  of  the  calibration  curve  [V/K], 


Figure  39  shows  the  deviation  as  a  function  of  AT.  The  deviation  is  less  1°C  even  for  AT  of 
100°C.  The  small  change  in  resistance  with  the  temperature  (Figure  37)  is  attributed  to  this 
behavior.  The  results  imply  that  the  junction  temperature  determined  from  the  forward  voltage  of 
the  LD  bar  at  the  operating  current  can  be  considered  as  the  average  junction  temperature  of  the 
LD  bar  in  practice.  This  implication  will  be  confirmed  later  with  the  actual  non-linear 
temperature  distribution  of  the  LD  bar. 
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Figure  39:  Deviation  of  the  Measured  Junction  Temperature  of  LD  Bar  from  the  Average 
Junction  Temperature  of  LD  Bar  with  a  Linearly  Changing  Temperature  of  AT  from  the 

Edge  to  the  Center 

3.2.3.3  Average  Junction  Temperature  Measurement 

Figure  40  shows  the  transient  voltage  behavior  of  the  LD  bar  obtained  after  blocking  the 
operating  current  of  80  A.  An  extreme  voltage  peak  at  the  beginning  of  transient  behavior  is 
clearly  visible.  It  was  produced  by  an  inductor  voltage  attributed  to  the  large  rate  of  current  change 
(79.88  A  to  zero)  and  the  non-zero  inductance  of  the  LD  [25,  26],  The  peak  was  large  but 
disappeared  quickly  after  200  ps. 
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Time  (ms) 

Figure  40:  Transient  Voltage  Behavior  of  the  LD  Bar  Obtained  after  Blocking  the 

Operating  Current  of  80  A 

Based  on  theoretical  analysis  [27-30],  it  is  known  that  a  junction  temperature  changes  linearly  in 
the  square  root  of  the  time  scale,  if  heat  is  dissipated  in  one  direction  through  a  homogenous 
material.  In  the  case  of  this  epi-down  LD  bar  with  the  water-cooled  microchannel  cooler,  the 
heat  transfer  through  the  GaAs  substrate,  and  convection,  as  well  as  radiation,  to  the  ambient 
surroundings  is  practically  negligible  (less  than  1%  of  the  total)  due  to  the  extremely  large  heat 
transfer  coefficient  in  the  water-cooled  microchannel  cooler  (this  will  be  discussed  further  in 
Section  3.2.4.2).  Thus,  the  linear  extrapolation  in  the  square  root  time  scale  is  applicable  for  the 
LD  bar.  The  enlarged  view  of  the  region  marked  by  a  dashed  box  in  Figure  41  is  shown  in 
Figure  41a.  The  voltage  was  converted  into  the  temperature  using  the  calibration  curve  and  it 
was  plotted  in  the  square  root  time  scale  (Figure  41b). 
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Time1c  {ms1/2} 

(b) 

Figure  41:  (a)  Enlarged  View  of  the  Region  Marked  by  a  Dashed  Box  in  Figure  6  and  (b) 
Average  Junction  Temperature  in  the  Square  Root  Time  Scale 

The  linear  extrapolation  provides  the  estimated  average  junction  temperature  at  the  operating 

current. 


The  transient  junction  temperature  behavior  can  be  divided  into  three  zones.  Zone  1  is  the 
region  dominated  by  the  electrical  delay.  Zone  2  is  the  region  where  the  linear  junction 
temperature  variation  follows  the  square  root  time  scale.  When  the  propagating  thermal  wave 
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reaches  the  microcooler  interface,  the  transient  junction  temperature  behavior  of  the  CuW 
submount  vanishes  and  we  enter  Zone  3.  It  is  estimated  from  Figure  7b  that  Zone  2  ends  at 
t  =1.57  ms  (=  1.25  ms1/2).  The  following  analytical  analysis  was  conducted  to  confirm  Zone  2. 

The  transient  domain  governed  by  the  CuW  submount  can  be  calculated  analytically  using  a  time 
constant,  which  can  be  expressed  as  [31,  32]: 

pc  d- 
r*  = — 1 — 


where  m  is  the  thermal  time  constant  [s];  d  is  the  thickness  along  the  heat  transfer  direction  [m]; 
k  is  the  thermal  conductivity  [W/m-K];  cp  is  the  specific  heat  [J/kg-K];  and  p  is  the  density 
[kg/m3].  Material  properties,  thickness,  and  the  calculated  thermal  time  constant  of  CuW  are 
listed  in  Table  2.  A  thermal  time  constant  value  for  the  CuW  was  determined  as  1.37  ms,  which 
is  defined  as  the  heating  or  cooling  time  required  to  produce  a  temperature  change  at  the  heat 
source  (junction)  equal  to  63.2%  of  the  total  temperature  difference  between  the  initial  and  the 
final  temperature.  This  value  is  reasonably  close  to  the  experimental  observation,  which  confirms 
the  validity  of  the  experimental  data. 


Table  2.  Material  Properties,  Thickness,  and  Calculated  Time  Constant  used  in  the 

Analytical  Solution  [19, 20] 


Material 

Density 

(kg/m3) 

Specific  heat 
(J/kg-K) 

Conductivity 

(W/m-K) 

Thickness 

(pm) 

Time  constant 
(ms) 

CuW 

17300 

160 

200 

300 

1.37 

The  average  junction  temperature  at  the  operating  current  was  estimated  from  the  linear 
extrapolation  shown  in  Figure  7b;  the  estimated  average  junction  temperature  at  80  A  was  35.6°C. 
The  discrepancy  in  the  repeatability  of  the  average  junction  temperature  measurement  was  less 
than  0.1  °C,  which  was  attributed  to  the  probe  current  source  inaccuracy. 

The  average  junction  temperatures  were  measured  from  10  A  to  80  A  at  an  interval  of  10  A.  The 
results  are  shown  in  Figure  42.  The  connected  lines  between  the  measured  data  represent  the 
trend  of  the  results.  As  expected,  the  junction  temperature  increased  with  the  current,  but  the  rate 
started  to  decrease  around  the  threshold  current  (26  A),  where  the  stimulated  emission  began  to 
occur;  i.e.,  the  higher  wall-plug  efficiency  lead  to  the  reduction  of  the  heat  dissipation  as  well  as 
the  junction  temperature.  The  heat  dissipation  as  a  function  of  the  forward  current  will  be 
discussed  further  in  Section  3.2.4. 1. 
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Figure  42:  Average  Junction  Temperature  at  different  Forward  Currents 

3.2.4  Heat  Dissipation  and  Microcooler  Effective  Heat  Transfer  Coefficient 

The  forward  voltage  and  the  emitted  radiant  flux  are  measured  to  quantify  the  amount  of  heat 
dissipation,  using  the  following  relationship  [33]: 

Ph(T1)  =  IrVf(T1)-<KTj)  (6) 

where  Ph  is  the  heat  dissipation  [W];  If  is  the  forward  current  [A];  Vf  is  the  forward  voltage  [V]; 
and  ®  is  the  radiant  flux  [W].  After  measuring  the  heat  dissipation  in  Section  5.1,  the  effective 
heat  transfer  coefficient  of  the  water-cooled  microchannel  is  calculated  inversely  from  the 
numerical  simulation,  using  the  measured  average  junction  temperature  and  the  heat  dissipation 
in  the  next  section. 

3.2.4. 1  Measurement  of  Heat  Dissipation 

The  DAQ1  measured  the  forward  voltage  of  the  LD  bar  and  the  optical  power  sensor  measured 
the  optical  power  (the  radiant  flux)  continuously.  When  the  optical  power  and  the  forward 
voltage  of  the  LD  reached  the  steady  state  condition,  the  values  of  forward  voltage  and  optical 
power  were  recorded,  from  which  the  heat  dissipation  was  calculated  using  Equation  (6). 

The  forward  voltages,  electrical  input  power  (product  of  the  forward  voltage  and  the  forward 
current),  optical  powers,  and  heat  dissipations  were  measured  as  a  function  of  current  with  an 
interval  of  10  A.  The  results  are  shown  in  Figure  43.  The  forward  voltage  and  the  electrical  input 
power  increased  with  the  current.  The  optical  power  was  virtually  negligible  before  the  threshold 
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current  (26  A)  and  increased  linearly  with  the  operating  current  after  the  threshold  current. 
Similar  to  the  junction  temperature,  the  heat  dissipation  increased  with  the  current  and  the  rate 
started  to  decrease  around  the  threshold  current. 


Figure  43:  (a)  Forward  Voltage  and  (b)  Electrical  Input  Power,  Radiant  Flux,  and  Heat 
Dissipation  as  a  Function  of  Forward  Current 
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The  thermal  resistances  can  be  estimated  by  dividing  the  temperature  difference  between  the 
average  junction  temperature  ( Tjave )  and  the  inlet  water  temperature  ( Tiniet=  20°C)  by  the  heat 
dissipation.  The  average  junction  temperatures,  the  heat  dissipations,  and  the  calculated  thermal 
resistances  (Rth)  are  summarized  in  Table  3.  The  results  showed  a  consistent  thermal  resistance 
of  0.295  K/W  ±0.015  K/W  for  all  measured  currents.  The  measurement  uncertainty  of  the  heat 
dissipation  was  mainly  caused  by  the  accuracy  of  the  optical  power  sensor. 


Table  3.  Average  Junction  Temperature  and  Thermal  Resistance  Estimations  at  different 
Forward  Currents  at  20°C  of  the  Inlet  Water  Temperature 


//[A] 

10 

20 

30 

40 

50 

60 

70 

80 

P*[  W] 

13.4  ± 

0.3 

27.4  ± 
0.4 

37.0  ± 

0.7 

39.2  ± 

1.2 

41.7  ± 

1.6 

44.5  ± 
2.2 

47.7  + 

2.7 

51.1  + 

3.1 

rr  pci 

23.8 

27.6 

30.2 

30.9 

31.6 

32.7 

34.1 

35.6 

ir  - 

PC] 

3.8 

7.6 

L0.2 

10.9 

11.6 

12.7 

14.1 

15.6 

if,*  [K/W] 

0.28  ± 
0.01 

0.28  ± 
001 

0.28  ± 
0.01 

0.28  ± 

0.01 

0.28  ± 
0.01 

0.29  ± 
0.01 

0.30  + 
0.02 

0.31  + 
0.02 

3.2.4.2  Effective  Heat  Transfer  Coefficient 

The  numerical  model  (ANSYS  16.1)  used  in  the  current  analysis  is  shown  in  Figure  44.  The 
model  has  the  same  geometry  as  the  LD  bar  (23  emitters;  emitter  width  of  200  pm  and  the  fill- 
factor  of  50%,  i.e.,  the  pitch  of  400  pm).  A  metallization  layer  (Ti/Pt/Au)  layer  between  the 
emitter  and  the  AuSn  solder  was  not  considered  in  the  model  due  to  the  ignorable  thermal 
resistance.  The  values  of  thermal  conductivity  of  GaAs,  CuW  submount,  and  AuSn  solder  used 
in  the  analysis  are  54  W/m-K  [34],  200  W/m-K  [19,  20],  and  58  W/m-K  [35],  respectively. 


Figure  44:  3D  Model 


The  ambient  temperature  was  set  at  20°C  (the  same  as  the  inlet  water  temperature).  The 
effective  heat  transfer  coefficient  of  the  water-cooled  microchannel  cooler  was  assumed  to  be 
uniform  on  the  bottom  of  the  CuW  submount.  The  effective  heat  transfer  coefficients  for 
natural  convection  (5  W/m2-K)  and  radiation  (GaAs  emissivity  of  0.62  [36])  were  set  on  the  top 
and  the  sides  of  the  model,  albeit  with  the  expectation  of  negligible  effects  on  the  junction 
temperature. 
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It  is  important  to  note  that  Equation  (6)  is  applicable  only  when  the  junction  temperature  is 
uniform.  In  addition,  the  uniform  heat  dissipation  would  be  desired  to  determine  the  effective 
heat  transfer  coefficient  most  accurately.  Thus,  the  lowest  operating  current  (10  A)  was  used  to 
calculate  the  effective  heat  transfer  coefficient. 

The  heat  dissipation  obtained  in  section  3.2.4. 1  was  applied  uniformly  on  the  emitters,  and  then, 
the  effective  heat  transfer  coefficient  was  adjusted  until  the  difference  between  the  measured 
average  junction  temperature  and  the  numerically  calculated  average  junction  temperature 
reached  its  minimum  value.  The  average  junction  temperature  difference,  after  typically  5 
iterations,  was  less  than  0.1  °C,  and  the  resulting  effective  heat  transfer  coefficient  was  found  to 
equal  98  kW/m2K. 

The  junction  temperature  distribution  at  20  A  was  also  calculated  to  validate  the  effective  heat 
dissipation.  The  difference  between  the  average  junction  temperatures  (experimental  and 
numerical)  was  less  than  0.1  °C,  which  confirmed  the  validity  of  the  effective  heat  dissipation. 

3.2.5  Numerical  Prediction  of  Junction  Temperature  Distribution 

It  is  important  to  understand  the  effect  of  the  junction  temperature  on  the  heat  dissipation  before 
performing  numerical  analyses  at  high  operating  currents  because  Equation  (6)  is  applicable  only 
when  the  junction  temperature  is  uniform.  The  heat  dissipation  at  80  A  was  measured  with  three 
inlet  water  temperatures  (10,  15,  and  20°C).  The  junction  temperature  at  each  inlet  temperature 
was  also  measured,  and  the  results  are  summarized  in  Table  4. 

Table  4.  Heat  Dissipation  at  80A  under  different  Inlet  Water  Temperatures 


r°c] 

10 

15 

20 

rr  ra 

25.6 

29.2 

35.6 

VfWl 

1.458 

1.457 

1.456 

I/Vf[  W] 

116.66 

116.58 

116.46 

<t>rwi 

67.80 

66.85 

65.34 

_ mu _ 

48.9 

49.7 

51.1 

As  the  average  junction  temperature  increased  from  25.6  to  35.6°C,  the  forward  voltage  as  well 
as  the  radiant  flux  decreased.  The  forward  voltage  reduction  reduced  the  total  electrical  power 
consumption,  while  the  radiant  flux  reduction  increased  the  fraction  of  the  input  power  converted 
to  heat.  It  is  worth  noting  that  the  net  heat  dissipation  increased  only  by  2.2  W  (or  4%) 
corresponding  to  the  junction  temperature  increase  of  only  0.6°C,  as  the  two  parameters 
compensated  their  effects  on  heat  dissipation  [37],  The  results  indicate  that  the  junction 
temperature  dependency  on  the  heat  generation,  over  a  range  of  10°C,  is  not  significant,  which 
provides  a  technical  rationale  for  the  following  numerical  study. 

3.2.5. 1  Temperature  Distribution  in  LD  Bar 

The  heat  dissipation  and  effective  heat  transfer  coefficient  obtained  in  the  previous  section  were 
used  to  predict  the  junction  temperature  distribution  over  the  LD  bar  numerically  under  high 
operating  currents.  The  junction  temperature  distribution  of  the  LD  bar  at  80  A  is  shown  in 

49 

Approved  for  public  release;  distribution  is  unlimited. 


Figure  45a.  The  GaAs  substrate  was  not  shown  in  order  to  clearly  show  the  temperature 
distribution  of  the  emitters. 

The  highest  temperature  occurred  at  the  emitting  side  of  the  center  emitter  and  the  lowest 
temperature  was  observed  at  the  opposite  side  of  the  edge  emitter.  The  junction  temperature 
decrease  toward  the  edge  emitter  and  the  back  end  was  attributed  to  the  effect  of  the  hear 
spreader.  Because  the  edge  emitter  and  the  back  end  had  cooling  area  enhancement  from  the 
CuW  submount,  the  more  heat  could  be  dissipated  due  to  the  extra  heat  spreading.  The 
maximum  temperature  difference  was  10.1°C.  This  result  confirms  that  the  application  of 
uniform  heat  dissipation  across  all  the  emitters  in  a  multi-emitter  LD  bar  can  be  expected  to 
provide  acceptable  numerical  results  for  input  current  of  up  to  80  A. 

The  junction  temperature  distribution  was  also  calculated  at  the  operating  current  of  160  A. 

Since  the  test  apparatus  was  only  capable  of  providing  125  A,  an  additional  power  supply 
(N5744A:  Keysight  Technologies)  was  connected  to  the  operating  current  source  in  parallel  to 
provide  the  additional  current  of  35  A.  The  forward  voltage  and  the  optical  power  at  160  A  were 
1.544  V  and  162.4  W,  respectively,  and  the  heat  dissipation  was  84.7  W  (65.7%  of  the  wall-plug 
efficiency).  It  is  to  be  noted  that  the  average  junction  temperature  could  not  be  measured  at 
160  A  because  the  threshold  current  of  the  MOSFET  switch  was  around  90  A. 

The  junction  temperature  distribution  of  the  LD  bar  at  160  A  is  shown  in  Figure  45b.  The  highest 
and  lowest  temperatures  are  47.9°C  and  31.7°C;  the  junction  temperature  difference  is  16.2°C. 


Figure  45:  Temperature  Distribution  of  the  LD  Bar  at  (a)  80  A  and  (b)  160  A 
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The  average  junction  temperature  of  each  emitter  is  compared  in  Figure  46,  where  the  only  left 
half  is  shown  due  to  the  symmetry.  The  average  junction  temperatures  remain  virtually 
unchanged  over  the  center  half  of  emitters  (from  #12  to  #7)  and  then  rapidly  drops  toward  the 
edge  emitter.  The  maximum  average  junction  temperature  differences  among  the  12  emitters  are 
4.9°C  and  7.8°C  at  80  A  and  160  A,  respectively. 


Figure  46:  Average  Junction  Temperature  of  each  Emitter  in  Left  Half  (symmetry) 

The  front- to-back  junction  temperature  variations  within  the  center  (#12)  and  the  edge  (#1) 
emitter  are  plotted  in  Figure  47.  The  junction  temperature  decreases  exponentially  from  the 
emitting  side  to  the  back  end.  The  junction  temperature  variations  within  the  emitter  are  largest 
at  the  center  emitter,  and  the  magnitude  increases  as  the  operating  current  increases.  The  junction 
temperature  variations  of  the  center  emitter  are  5.4°C  and  8.7°C  for  80  A  and  160  A, 
respectively. 
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Figure  47:  Junction  Temperature  Variations  along  the  Emitter 

The  deviations  of  the  true  junction  temperature  from  the  average  junction  temperature  of  the  bar 
were  calculated  with  the  actual  non-linear  temperature  distribution  of  the  LD  bar.  Based 
Equation  3  was  first  used  to  determine  the  true  forward  voltage  and  the  forward  voltage  from  the 
average  junction  temperature: 


V*™  (80)  =  1 . 1 682  3  9  V  and  V™ (80)  =  1 , 1 68242  V .  md 
V*'“e (1 60)  =  1 . 1 57465  V  and  V™  (160)  =  1.15  7472  V 


The  deviations  were  then  calculated  from  Equation  4;  they  were  0.002°C  and  0.006°C  for  80  A 
and  160  A,  respectively.  As  expected  from  the  small  change  in  resistance  with  the  temperature, 
the  deviations  were  negligible.  The  results  confirmed  the  validity  of  the  proposed  method. 

3.2.5.2  Wall-plug  Efficiency  and  Spectral  Power  Distribution 

The  junction  temperature  change  affects  the  wall-plug  efficiency,  which  is  defined  as  the  optical 
output  power  divided  by  the  electrical  input  power.  Each  emitter  of  the  LD  bar  tested  in  the 
study  produced  an  optical  power  of  »7.06  W  at  the  operating  current  of  6.95  A  (the  LD  bar 
with  23  emitters  at  the  operating  current  of  160  A). 

It  was  reported  in  Ref.  [38]  that  for  a  single  emitter  (a  center  wavelength  of  975  nm) 
producing  6  W  at  6A,  the  wall-plug  efficiency  increased  by  only  3%  from  275  K  (71%)  to  260  K 
(74%).  Thus,  it  is  reasonable  to  assume  that  each  emitter  of  the  LD  bar  tested  in  the  study  has 
similar  behavior  to  that  of  the  single  emitter  reported  in  Ref.  [38],  and  consequently,  the  effect  of 
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the  junction  temperature  distribution  on  the  wall-plug  efficiency  is  not  significant. 


The  power  spectrum  of  the  LD  bar  was  obtained  by  a  spectrometer  (AvaSpec-ULS3648-  USB2: 
Avantes)  combined  with  a  cosine  corrector  (CC-UV/VIS/NIR-5.0:  Avantes).  The  diameter  of  the 
corrector  was  larger  than  the  LD  bar  to  ensure  that  the  spectrometer  received  light  from  all  23 
emitters  equally.  The  output  obtained  from  the  measurement  system  was  the  irradiance  (W/m2) 
of  the  LD  bar  as  a  function  of  wavelength.  The  power  spectrum  was  obtained  by  normalizing  the 
irradiance  distribution  by  the  peak  irradiance.  The  results  obtained  at  30,  60,  and  120  A  are 
shown  in  Figure  48. 


>  I  r  i  i — T~  T  |  i  >-  i  | 

924  926  928  930  932  934 


Wavelength  (nm) 

Figure  48:  Normalized  Power  Spectrum  at  30,  60,  and  160  A 


LDs  are  known  to  show  the  spectral  red  shift  (i.e.,  higher  peak  wavelength)  at  higher  junction 
temperatures  due  to  the  reduced  band  gap  energy  [39].  The  spectral  red  shift  caused  by  the 
junction  temperature  has  been  reported  to  be  0.28  nm/K  [40]  and  0.32  nm/K  [41]  for  808  nm  and 
980  nm  LDs,  respectively.  The  results  in  Figure  14  show  the  red  shift  of  «0.3  nm/K,  which  is 
consistent  with  the  reported  values. 


The  results  also  show  that  the  spectrum  is  broadened  as  the  current  increases;  i.e.,  the  larger  full 
width  at  half  maximum  (FWHM)  at  the  higher  forward  current.  The  asymmetry  of  the  spectrum 
(i.e.,  more  broadening  toward  the  lower  wavelength)  become  more  severe  as  the  current 
increases.  Both  broadening  and  asymmetry  are  attributed  to  the  larger  junction  temperature 
gradient  at  the  high  currents. 

The  analysis  was  based  on  a  very  large  heat  transfer  coefficient.  In  practice,  various  thermal 
solutions  can  be  employed  for  cooling  the  LD  bar.  In  terms  of  the  coefficient  of  performance 
(COP),  the  lower  heat  transfer  coefficients  (i.e.,  the  decreased  flow  rate)  can  reduce  the  operating 
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costs.  However,  the  junction  temperature  will  rise  and  the  temperature  variation  within  the  LD 
bar  will  also  increase  with  lower  heat  transfer  coefficients.  This  will  increase  the  asymmetry  of 
the  SPD  and  the  peak  wavelength  shift,  which  can  reduce  the  pumping  efficiency  [42]. 

A  high  heat  transfer  coefficient  is  desired  to  increase  the  pumping  efficiency,  which  can  be 
achieved  with  an  extreme  flow  rate  of  a  coolant.  However,  the  higher  flow  rate  reduces  the  COP, 
which  increases  the  operating  cost.  In  addition,  it  can  accelerate  the  erosion  process  of  the 
surface  structures  inside  microchannels,  which  will  increase  the  junction  temperature  and  will 
eventually  reduce  the  lifetime  of  the  LD  bar  [43,  44].  Consequently,  optimization  of  thermal 
solutions  for  high  power  LD  bars  should  be  sought  while  considering  the  operating  cost  as  well 
as  various  thermal,  mechanical,  and  optical  aspects  of  the  system. 

3.2.6  Summary 

A  hybrid  experimental  and  numerical  method  was  proposed  and  implemented  for  predicting  the 
junction  temperature  distribution  of  a  high  power  LD  bar.  A  commercial  water-cooled  LD  bar 
was  utilized  to  illustrate  and  validate  the  proposed  method.  The  average  junction  temperature  and 
the  heat  dissipation  were  measured,  and  the  effective  heat  transfer  coefficient  of  the  cooling 
system  was  determined  inversely  using  numerical  simulation.  The  characterized  properties  were 
used  to  predict  the  junction  temperature  distributions  of  the  LD  bar  at  the  extreme  operating 
currents.  The  results  showed  significant  junction  temperature  variations  not  only  among  emitters 
(7.8°C)  but  also  along  each  emitter  (8.7°C)  at  160A,  which  increased  the  asymmetry  of  the 
power  spectrum.  The  proposed  method  can  be  used  to  determine  the  proper  operation  condition 
of  the  LD  bar  as  well  as  to  evaluate  designs  during  packaging  platform  development.  The  future 
work  will  address  a  methodology  to  define  the  optimum  design  of  LD  bars  considering  the  COP, 
performance,  and  reliability. 

3.3  Spectral  Power  Distribution  Deconvolution  Scheme  for  High-Power  LD  Bar 
3.3.5  Introduction 

A  typical  LD  array  of  10  mm  length  contains  10-60  emitters  with  fill  factors  from  10%  to  90%. 
The  lateral  heat  spreading  in  an  LD  array  causes  the  thermal  crosstalk  effect  between  emitters, 
and  thus,  the  junction  temperature  distribution  typically  has  a  large  variation  [51].  In  addition, 
hotter  emitters  in  an  LD  array  take  a  larger  share  of  the  total  array  current,  and  emit  more  optical 
power.  The  effect  is  known  as  “current  competition”,  which  causes  the  non-uniform  power 
distribution  [45,  46],  Consequently,  the  SPDs  of  individual  emitters  would  have  significant 
variations. 

The  SPDs  of  an  LD  array  with  non-uniform  temperature  and  power  distributions  are  illustrated  in 
Figure  49,  where  SPDAnay,  SPD1  and  SPD4  represent  the  SPDs  of  the  LD  array,  the  edge  emitter 
and  the  center  emitter,  respectively.  The  junction  temperature  of  the  center  emitter  is  always 
higher  compared  to  the  edge  emitter.  As  a  result,  the  SPD  of  the  center  emitter  (SPD4)  would 
have  the  maximum  power  and  the  longest  wavelength,  while  the  SPD  of  the  edge  emitter  (SPD1) 
the  lowest  power  and  the  shortest  wavelength. 
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The  SPD  of  an  LD  array  can  be  readily  measured  by  a  spectrometer  connected  with  a  cosine 
corrector.  The  SPDs  of  individual  emitters  can  also  be  measured  by  placing  a  beam  baffle  in 
front  of  the  array.  In  practice,  however,  the  distance  between  the  beam  baffle  and  the  LD  array 
has  to  be  extremely  small  because  of  the  small  pitch  between  adjacent  emitters  and  the  large 
beam  divergence  of  the  LD  array.  The  optical  feedback  from  the  beam  baffle  at  such  a  close 
distance  can  cause  degradation  or  even  catastrophic  optical  mirror  damage  of  the  emitters  [9,10], 
In  addition,  translating  the  baffle  accurately  to  open  only  one  emitter  while  keeping  the  baffle  at 
the  close  proximity  of  the  array  is  very  challenging  in  practice,  especially  for  LD  arrays  with 

high  fill  factors  (>50%). 


k 


Figure  49:  Illustration  of  SPDs  for  an  LD  Array  with  a  Non-Uniform  Temperature  and 

Power  Distributions 

Where  SPDArray,  SPD1,  and  SPD4  represent  the  SPDs  of  an  LD  array,  the  edge  emitter  and  the 

center  emitter,  respectively 

The  objective  of  this  paper  is  to  propose  a  novel  method  to  predict  the  SPDs  of  individual 
emitters  by  deconvoluting  the  SPD  of  an  LD  array.  The  proposed  method  takes  into  account  the 
thermal  cross  talk  effect  as  well  as  the  current  competition  effect.  A  complete  analytical 
description  of  the  proposed  method  is  described  in  Section  3.  The  implementation  of  the 
proposed  method  using  a  commercial  LD  array  is  presented  in  Section  4.  Applications  of  the 
SPD  deconvolution  are  presented  in  Section  5. 

3.3.2  Analytical  Model  for  Deconvolution 

The  SPD  of  an  LD  array,  SPDARRAY(AJ,  is  simply  the  sum  of  the  SPDs  of  single  emitters,  and  can 
be  expressed  as: 

N 

.SPD  ^  (A)  =  V  SPD1  (A) 
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where  SPD1  (A),  is  the  SPD  of  the  ith  emitter;  and  N  is  the  number  of  emitters. 

Figure  50  illustrates  schematically  the  SPD  of  a  single-emitter  in  an  LD  array.  The  maximum 
power  and  the  full  FWHM  of  the  ith  emitter  are  denoted  as  PW  and  W,  and  p  is  the  central 
wavelength  of  the  zth  emitter,  which  is  defined  as  the  wavelength  satisfying  thet  condition  that 
is  equal  to  at  ^ax/2. 


Figure  50:  Schematic  Illustration  of  the  SPD  of  a  Single-Emitter  in  an  LD  Array 

In  a  typical  LD  array,  an  emitter  located  at  the  center  of  an  LD  array  (referred  to  as  “center 
emitter”)  has  the  highest  junction  temperature.  Its  normalized  SPD  can  be  defined  as: 


-center 


SPD  (A)  = 


SPD""*' (A) 


Q/cmter 

max 


(8) 


where  SPDcenter  (A)  is  the  SPD  of  the  center  emitter  normalized  by  the  maximum  power  of  the 
center  emitter  spectrum,  . 

The  SPD  of  the  center  emitter  is  asymmetric  because  the  profile  of  the  gain  spectrum  is  not 
symmetric  [47].  The  asymmetric  SPD  can  be  expressed  using  multiple  Gaussian  functions: [48, 
49] 
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where  k  is  the  number  of  Gaussian  functions;  Aj  is  the  normalized  amplitude  of  the  /h  Gaussian 
function;  A,  is  the  central  wavelength  of  the  /h  Gaussian  function;  and  w/  is  proportional  to  the 
FWHM  (Wj)  of  the /h  Gaussian  function,  which  is  defined  as: 

...  5 

J  2An2 

The  junction  temperature  of  each  emitter  varies  within  an  LD  array.  It  has  been  known  that  the 
central  wavelength  of  the  SPD  changes  linearly  with  the  temperature. [47]  As  the  junction 
temperature  increases,  the  band-gap  energy  decreases,  and  the  refractive  index  and  cavity  length 
increase.  The  band-gap  energy  reduction  is  the  dominant  factor  causing  the  central  wavelength 
shift.  The  relationship  between  the  central  wavelength  shift  and  the  junction  temperature,  then, 
can  be  expressed  as: 


a  i  i  i  i  cent®'  {  T1  /  -reenter  \ 

A A  =  AC-AC  =a[T  -T  ) 


(ID 


i  center 

where  xc  and  A  are  the  central  wavelengths  of  the  i  emitter  and  the  center  emitter  T  and 
jcmter  are  junction  temperature  of  the  ith  emitter  and  the  center  emitter;  and  a  is  the 
temperature  coefficient  of  wavelength  that  can  be  determined  experimentally. 


As  described  in  Refs.  [45,  46],  the  maximum  power  of  each  emitter  changes  with  temperature  by 
the  effect  known  as  “current  competition”.  The  emitters  that  turn  on  earlier  take  a  larger  share  of 
the  total  array  current,  and  emit  more  power  because  hotter  emitters  in  an  LD  array  have  a 
reduced  bandgap  energy,  and  thus  a  lower  threshold  current. 


Based  on  the  theoretical  and  experimental  results  by  S.  Bull  et  al.  [45,  46],  the  effect  of  current 
competition  can  be  approximated  described  by  an  exponential  function.  The  maximum  powers  of 
the  ith  emitter  and  the  center  emitter,  then,  can  be  expressed  as: 


P1 

max 


(12) 


^center  _ ^ 

max  ^ 


(13) 


where  B  and  C  are  constants.  By  combining  Equations  (12)  and  (13),  the  maximum  power  of  the 
ith  emitter  normalized  by  can  be  described  as: 
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p* 

max 

■pcentsr 

^ItHX 


(14) 


The  constant,  B,  will  be  referred  to  as  the  current  competition  constant  that  can  be  determined 
experimentally. 

It  was  reported  that  the  shape  of  the  normalized  gain  spectrum  profile  remains  virtually  the  same 
regardless  of  junction  temperatures  [50],  This  implies  that  the  FWHM  of  the  normalized  single 
emitter  SPD  will  not  be  altered  by  the  junction  temperature  variations  with  an  LD  array.  Using 
Equations  (9)  and  (1 1),  the  SPD  of  the  ith  emitter,  then,  can  be  expressed  as: 


SPD'iX)  =  ■  SPD  ^ (A- AX)  = 
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Using  Equation  (15),  Equation  (16)  can  be  written  as: 
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3.3.3  Implementation 

In  order  to  predict  the  SPDs  of  individual  emitters  of  an  LD  array  using  Equitation  (16),  four 
key  parameters  have  to  be  determined  experimentally:  (1)  the  normalized  SPD  of  the  center 
emitter,  (2)  the  temperature  coefficient  of  wavelength,  a,  (3)  the  current  competition  constant,  B, 
and  (4)  the  maximum  power  of  the  center  emitter,  F„‘X'.  The  junction  temperature  distribution  of 
the  array  is  critically  required  to  determine  the  four  parameters.  The  hybrid  experimental/numerical 
method  proposed  previously  by  the  authors  [51]  is  employed  to  determine  the  required  junction 
temperature  distribution.  Detailed  procedures  to  determine  the  four  parameters  are  presented  in 
sections  333.3  to  3. 3. 3.6,  respectively,  after  describing  the  testing  apparatus  and  the  hybrid 
method  for  the  junction  temperature  measurement  in  sections  3.3.3. 1  and  3. 3. 3.2. 
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3.3.3. 1  Device  and  Measurement  Apparatus 


A  commercial  930  nm  LD  array  (E11.4N-940.10-150C-S013.1:  DILAS)  used  in  the  study 
consists  of  23  identical  emitters.  The  width  of  each  emitter  is  200  pm  and  the  pitch  between 
adjacent  emitters  is  400  pm  (a  fill-factor  of  50%).  The  maximum  optical  power  at  160  A  is 
160  W.  The  LD  chip  is  epi-down  bonded  on  a  CuW  submount  using  AuSn  die  attach.  The  CuW 
submount  (CTE):  6.5  ppm/°C)  is  placed  between  the  GaAs  chip  (6.4  ppm/°C)  and  the  copper 
microchannel  heat  sink  (16.6  ppm/°C)  to  reduce  the  mismatch  in  the  CTE  [18,  52], 


A  test  apparatus  to  measure  the  SPD  of  the  central  emitter  is  illustrated  in  Figure  5 1 .  The  power 
supply  (LDX-36125-12:  ILX  Lightwave)  applies  the  operating  current  with  a  nominal  accuracy 
of  ±0.1%.  The  spectrometer  (AvaSpec-3648)  connected  with  a  cosine  corrector  measures  the 
spectrum.  The  wavelength  range  of  the  spectrometer  is  from  200  nm  to  1100  nm,  and  the 
resolution  is  0.025  nm.  The  chiller  (ISOTEMP  1 115V/60HZ  PD-1:  Fisher  Scientific)  regulates 
the  inlet  water  temperature  with  a  temperature  stability  of  ±0.1  °C.  The  flow  meter 
(FLDW321 1G:  OMEGA  Engineering)  controls  the  flow  rate  from  0  to  500  mL/min. 


LabVIEW 

controlled 

computer 


Spectrometer 


Top  view 
of  LD  bar 


Power  supply 
% 


Chiller 


Flow 

regulator 


Cosine 

corrector 


Figure  51:  Schematic  Illustration  of  the  SPD  of  a  Single-Emitter  in  an  LD  Array 


When  only  the  SPD  of  the  center  emitter  is  to  be  measured,  the  beam  baffle  is  placed  between 
the  LD  array  and  the  cosine  corrector.  The  baffle  is  made  of  graphite  to  minimize  the  optical 
feedback  while  dissipating  the  heat  generated  by  the  light  beam  effectively. 


In  order  to  measure  the  absolute  optical  power,  the  cosine  corrector  and  the  spectrometer  are 
replaced  to  the  optical  power  sensor  (USB-PM- 150-50:  Coherent  Laser  Group).  The  power 
supply,  the  DAQ  (USB-6212:  National  Instruments),  and  the  optical  power  sensor  are  integrated 
into  a  LABVIEW  program. 


3.3.3.2  Determination  of  Junction  Temperature  Distribution 

The  hybrid  experimental/numerical  method  [51]  was  developed  to  determine  the  junction 
temperature  distribution  within  a  high  power  LD  array.  With  the  method,  the  forward  voltage 
method  is  first  implemented  in  a  unique  experimental  setup  to  measure  the  average  junction 
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temperatures  of  the  LD  array.  After  measuring  the  heat  dissipation  of  the  LD  array,  the  effective 
heat  transfer  coefficients  of  the  cooling  system  at  different  flow  rates  are  determined  inversely 
from  the  numerical  simulation  using  the  measured  average  junction  temperature  and  the  heat 
dissipation.  The  characterized  effective  heat  transfer  coefficients  at  different  flow  rates  are  used 
to  predict  the  junction  temperature  distribution  over  the  LD  array  at  different  inlet  water 
temperatures.  More  details  about  the  method  can  be  found  in  the  reference  [51]. 

The  numerical  model  (ANSYS  Icepak  17.2)  used  in  the  analysis  is  shown  in  Figure  52.  The 
values  of  thermal  conductivity  of  GaAs,  CuW  submount,  and  AuSn  solder  are  54  W/m-K,  209 
W/m-K,  and  58  W/m-K,  respectively  [53-57].  The  ambient  temperature  was  set  to  be  20°C.  The 
effective  heat  transfer  coefficients  for  natural  convection  (5  W/(m2-K))  and  radiation  (GaAs 
emissivity  of  0.62)  were  set  on  the  top  and  the  sides  of  the  model,  although  they  had  negligible 
effects  on  the  junction  temperature  [58].  The  effective  heat  transfer  coefficient  of  the  water- 
cooled  microchannel  was  assumed  uniform  on  the  bottom  of  the  CuW  submount. 


Figure  52:  Schematic  Illustration  of  the  SPD  of  a  Single-Emitter  in  an  LD  Array 

The  generated  heat  from  the  active  region  of  the  LD  was  applied  uniformly  on  the  emitters  in  the 
LD  array,  and  the  effective  heat  transfer  coefficient  was  calculated  using  the  iterative  method. 

The  forward  voltage,  optical  power,  heat  dissipation,  average  junction  temperatures,  and 
calculated  effective  heat  transfer  coefficients  at  different  flow  rates  (100,  200,  300,  400,  and, 
500  mL/min)  are  summarized  in  Table  5.  Using  the  effective  heat  transfer  coefficients  and  heat 
dissipation,  the  temperature  distributions  of  the  LD  array  at  different  flow  rates  were  predicted. 
A  representative  steady-state  temperature  distribution  of  the  LD  array  is  shown  in  Figure  53, 
where //=  80  A,  h  =  80,500  W/(im-K),  and  Twet=  20°C.  The  maximum  temperature  occurred  at 
the  front  facet  of  the  center  emitter.  The  junction  temperature  decreased  towards  the  edge  emitter 
and  the  rear  facet  due  to  the  heat  spreading  effect  of  the  CuW  submount. 
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Table  5.  Effective  Heat  Transfer  Coefficients  at  different  Flow  Rates  with  an  Inlet  Water 

Temperature  of  20°C 


Flow  rate  [mL/min] 

100 

200 

300 

400 

500 

Forward  voltage  [V] 

1.4135 

1.4147 

1.4154 

1.4156 

1.4160 

Optical  power  [W] 

14.675 

16.215 

16.707 

16.986 

17.153 

Heat  dissipation  [W] 

42.034 

40.542 

40.079 

39.808 

39657 

Average  junction 
temperature  [°C] 

39.3 

34.9 

31.5 

31.0 

30.0 

Effective  heat  transfer 
coefficient  [kW/(m:  K)] 

56.25 

80.50 

112.65 

119.90 

137.60 

_  40.9  °C 
_ 

_  35.9  °C 
1 

30.8  °C 
■  25.7  °C 


Figure  53:  Temperature  Distribution  of  the  LD  Array 

(If  =  80  A,  h  =  80,500  W/(m2-K),  and  Twet=  20  °C) 


The  simulated  average  junction  temperature  of  each  emitter  at  different  heat  transfer  coefficients 
is  shown  in  Figure  54,  where  the  left  half  of  the  LD  array  is  shown  due  to  the  symmetry  of  the 
LD  array.  The  junction  temperatures  remain  nearly  unchanged  in  the  center  of  the  LD  array,  but 
rapidly  decrease  towards  the  edge  emitters.  As  expected,  the  temperature  variation  between  the 
center  emitter  and  the  edge  emitter  increases  with  the  low  heat  transfer  coefficient. 
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Figure  54:  Average  Junction  Temperature  of  each  Emitter  in  the  Left  Half  (symmetry) 

under  various  Heat  Transfer  Coefficients 

(If  =80  A,  andTiniet  =  20  °C). 


3.3.3.3  Determination  of  Normalized  SPD  of  Center  Emitter 

Considering  the  beam  divergence  along  slow  axis  (10°)  and  the  high  fill  factor  (50%)  of  the 
commercial  LD  array,  the  beam  baffle  width  and  the  distance  between  the  beam  baffle  and  the 
LD  array  should  be  smaller  than  0.36  mm  and  0.90  mm,  respectively  (Figure  55),  which  poses 
implementation  difficulties.  In  order  to  cope  with  the  problem,  a  larger  baffle  width  of  1.5  mm 
was  used  in  the  actual  experiments,  which  averaged  three  emitters  in  the  middle  of  the  array.  The 
larger  baffle  width  was  rationalized  by  the  fact  that  that  a  few  emitters  in  the  middle  of  the  array 
were  known  to  have  a  virtually  identical  junction  temperature. [51]  Two  supplementary 
experiments  were  conducted  to  confirm  the  rationale. 


Figure  55:  Illustration  of  the  Size  and  Location  of  an  Optical  Baffle  to  Measure  the  SPD  of 

a  Center  Emitter 
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In  the  first  experiment,  two  baffle  widths  of  1.5  and  3  mm  were  used  to  record  three  and  seven 
emitters  in  the  middle,  respectively.  The  normalized  SPDs  under  the  forward  current  (//)  of  80  A, 
the  flow  rate  if)  of  500  mL/min,  and  the  inlet  water  temperature  ( Tiniet )  of  20°C  are  compared  in 
Figure  56(a).  The  normalized  SPDs  are  nearly  identical,  which  clearly  indicates  that  even  seven 
emitters  in  the  center  have  the  virtually  same  SPDs  considering  the  measurement  uncertainty. 
Accordingly,  it  is  reasonable  to  conclude  that  the  SPD  obtained  from  three  center  emitters 
represents  the  normalized  SPD  of  a  single  emitter. 

In  the  second  experiment,  the  baffle  width  was  fixed  to  be  1 .5  mm  and  the  normalized  SPDs  of 
three  center  emitters  were  measured  at  different  cooling  conditions.  The  normalized  SPD  of  three 
center  emitters  at  two  conditions  are  compared  in  Figure  56(b):  flow  rates  and  inlet  water 
temperatures  of  (1)  300  mL/min  and  10°C  and  (2)  500  mL/min  and  20°C.  The  FWHMs  of  the 
SPDs  under  the  above  cooling  conditions  were  determined  from  the  SPDs.  Both  SPDs  have  the 
identical  FWHM  of  1.3  nm,  which  indicates  that  the  spectrum  width  of  the  center  emitters 
remains  unchanged  regardless  of  temperatures.  The  result  confirms  the  fact  [50]  that  the  shape  of 
the  normalized  gain  spectrum  profile  is  not  altered  by  the  junction  temperature  variation  within 
an  LD  array,  and  thus,  provides  a  technical  rationale  for  Equation  (15);  i.e.,  the  normalized  SPD 
of  each  emitter  in  an  LD  array  can  be  determined  from  the  normalized  SPD  of  the  center  emitter. 
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Wavelength  (nm) 


(b) 

Figure  56:  Normalized  SPDs  obtained  from  (a)  Three  and  Seven  Emitters  in  the  Middle  of 
an  LD  Array  where  the  Baffle  Widths  are  1.5  and  3  mm,  respectively  (//=  80  A ,/=  500 
mL/min,  and  Tiniet  =  20  °C)  and  (b)  Three  Emitters  in  the  Middle  of  an  LD  Array  at 
different  Flow  Rates  and  Inlet  Water  Temperatures 
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As  mentioned  earlier,  the  SPD  of  the  center  emitter  is  not  symmetric,  and  multiple  Gaussian 
functions  are  necessary  to  deconvolute  the  normalized  SPD  of  the  center  emitter.  [48]  Three 
Gaussian  functions  were  used  to  fit  the  normalized  SPD  of  the  center  emitter  in  this  study.  The 
results  are  shown  in  Figure  57,  where  the  measured  SPD  is  compared  with  the  Gaussian  fitting. 
Perfect  agreement  is  evident;  the  R2  value  is  close  to  1 .  The  Gaussian  fitting  parameters  are 
shown  in  Table  6. 


Figure  57:  SPD  Deconvolution  of  the  Normalized  SPD  of  the  Center  Emitter  using  Three 

Gaussian  Functions 

(If  =  80A,f=  200  mL/min,  and  Tiniet  =  20  °C) 

Table  6.  Three  Gaussian  Functions  used  to  define  the  Normalized  SPD  of  the  Center 

Emitter 


Ai  (nm) 

A2(nm) 

A3(nm) 

928.588 

929.114 

929.765 

wi  (nm) 

W2(nm) 

W3  (nm) 

0.613 

0.810 

0.293 

A\ 

a2 

A3 

0.468 

0.573 

0.177 
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3.3.3.4  Temperature  Coefficient  of  Wavelength 

The  SPDs  of  three  center  emitters  were  measured  at  various  junction  temperatures  to  determine 
the  temperature  coefficient  of  wavelength,  a.  The  normalized  SPDs  and  the  corresponding 
average  junction  temperatures  of  the  center  emitter  are  shown  in  Figure  58(a),  from  which  the 
central  wavelengths  can  be  obtained. 

The  central  wavelengths  of  center  emitter  are  plotted  as  a  function  of  average  junction 
temperature  in  Figure  58(b);  a  linear  relationship  is  evident.  The  temperature  coefficient  of 
wavelength  was  obtained  from  Figure  58(b)  using  Equation  (11);  it  was  0.3  nm/K.  The  spectral 
redshift  caused  by  the  junction  temperature  has  been  reported  as  0.26-0.28  nm/K  [59-61],  and 
0.32  nm/K  [62],  for  the  LDs  with  the  central  wavelength  of  808  nm  and  980  nm,  respectively. 
The  measured  value  is  consistent  with  the  reported  values. 
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(a) 


(b) 

Figure  58:  (a)  Measured  normalized  SPDs  representing  the  Center  Emitter  at  different 
Average  Junction  Temperatures  and  (b)  Central  Wavelengths  of  SPDs  plotted  as  a 
Function  of  the  Average  Junction  Temperature  of  the  Center  Emitters;  the  Linear 
Relationship  defines  the  Temperature  Coefficient  of  Wavelength 
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3.3.3.5  Determination  of  Current  Competition  Constant 

The  current  competition  constant,  B,  cannot  be  determined  deterministically  as  Equation  (14) 
contains  two  unknowns.  The  value  was  determined  iteratively.  An  initial  value  of  B  =  3.6  was 
first  estimated  by  fitting  the  data  in  Ref.  [46]  through  Equation  (14).  Then,  the  normalized 
SPD  of  each  emitter,  based  on  the  normalized  SPD  of  the  center  emitter  (assuming  P^*r=  1), 
was  calculated  using  Equation  (16),  and  subsequently  the  normalized  SPD  of  the  LD  array  was 
calculated  using  Equation  (7).  Finally,  a  non-linear  regression  was  performed  while  adjusting  the 
values  of  B  until  the  coefficient  of  determination,  R2,  between  the  simulated  result  and  the 
measured  normalized  SPD  of  the  LD  array  reached  its  maximum. 

The  normalized  SPD  of  the  LD  array  at  a  heat  transfer  coefficient  of  98,000  W/(m2  K)  (the 
operating  condition  recommended  by  the  LD  array  manufacturer)  with  an  inlet  water  temperature 
of  20°C  was  measured  to  determine  B.  The  results  from  the  iteration  process  are  shown  in 
Figure  59.  The  simulated  SPDs  of  each  emitter  obtained  from  the  initial  value  of  B  =  3.6  is  shown 
in  Figure  59(a).  The  predicted  SPDs  of  the  LD  array  with  the  initial  and  final  values  (B  =  3.6  and 
5.9)  are  compared  with  the  measured  normalized  SPDs  in  Figure  59(b)  and  (c),  respectively.  With 
the  correct  value  of  B,  the  predicted  SPD  became  virtually  identical  to  the  measured  SPD. 


(a) 
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Figure  59:  Predicted  (a)  Normalized  SPDs  of  Emitters  in  the  Left  Half  (symmetry),  (b) 
Normalized  SPD  of  the  LD  Array  using  an  Initial  Value  of  B  =  3.6,  and  (c)  Normalized  SPD 

of  LD  Array  using  the  Linal  Value  of  B  =  5.9 

(If  =  80  A,  h  =  98,000  W/(m*-K),  and  Tiniet  =  20  °C) 
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3.3.3.6  Determination  of  Maximum  Power  of  Center  Emitter 


From  Equations  (1)  and  (10),  the  SPD  of  the  LD  array  can  be  expressed  as: 


SPDim>(X)  =  Y^SPD'O.)  = 
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The  integration  of  SPD^rray  (A)  should  be  equal  to  the  total  optical  power  of  the  LD  array. 
Therefore,  the  maximum  power  of  the  center  emitter  P^**  ,  can  be  expressed  as: 
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where  Pa  is  the  total  optical  power  of  an  LD  array,  which  can  be  measured  experimentally. 

Using  the  measured  optical  power  of  the  LD  array,  the  values  of  P.PP"  under  various  cooling 
conditions  were  calculated  using  Equation  (18).  The  results  are  summarized  in  Table  7.  The 
value  of  PP^*"  decreases  with  the  increased  flow  rate,  which  is  attributed  to  more  uniform 
temperature  distributions  at  higher  flow  rates.  On  the  other  hand,  the  value  of  PP^‘  increases 
with  the  decreased  inlet  water  temperature  due  to  the  higher  optical  power  at  lower  junction 
temperatures. 


Table  7.  P^f  under  different  Cooling  Conditions 


Flowrate 

[nfL'min] 

500 

300 

100 

500 

300 

100 

500 

300 

100 

Inlet  water 
temperature 
PC] 

20 

20 

20 

15 

15 

15 

10 

10 

10 

c*  [w] 

2.21 

2.24 

2.40 

2,27 

2.29 

2.45 

2.30 

2.33 

2.47 
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3.3. 3. 7  Prediction  of  the  Absolute  SPDs 

Using  the  parameters  determined  in  the  previous  sections,  the  absolute  SPDs  of  all  emitters  were 
calculated  using  Equation  (16).  The  results  obtained  at  h  =  137,600  W/(m2-K)  and  Tiniet=  20°C 
are  shown  in  Figure  60(a).  The  seven  center  emitters  (#9,  #10,  #1 1  #12,  #13,  #14,  #15)  have  a 
similar  SPD  (i.e.,  the  maximum  amplitude,  the  central  wavelength,  and  the  spectrum  width),  as 
expected  from  the  experimental  results  reported  in  Section  3C.  The  SPD  of  the  edge  emitter  has 
a  lower  amplitude  and  a  shorter  wavelength  compared  with  the  center  emitter.  The  amplitude  of 
SPD  decreased  rapidly  towards  edge  emitter,  which  is  attributed  to  the  current  competition 
effect.  The  current  competition  effect  on  individual  emitters  will  be  discussed  further  in  the  next 
section. 

The  SPD  of  the  LD  array  at  h  =  137,600  W/(  m2K)  with  Tmiet  of  20°C  was  predicted  using 
Equation  (18).  The  result  is  compared  with  the  measured  SPD  in  Figure  60(b).  The  predicted 
SPD  agrees  well  with  the  experimental  data  in  both  shape  and  magnitude.  The  result 
corroborates  the  validity  of  the  proposed  method. 
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Wavelength  (nm) 

(b) 

Figure  60:  (a)  Predicted  Absolute  SPDs  of  Individual  Emitters  in  the  Left  Half  (symmetry) 
and  (b)  Comparison  between  Predicted  and  Measured  Absolute  SPD  of  the  LD  Array 

(If  =  80  A,  h  =  137,600  W/Um-K),  and  Tiniet  =  20  °C) 


3.3.4  SPDs  of  Single  Emitters  at  Different  Cooling  Conditions:  Results  and  Discussions 

The  SPDs  of  the  LD  array  obtained  at  various  cooling  conditions  are  shown  in  Figure  61.  They 
were  measured  at  hi  =  56,250,  hi  =  80,500,  hi  =  108,100,  h4  =  112,650,  and  hs  =  137,600 
W/(m2  K);  the  corresponding  flow  rates  were  100,  200,  300,  400,  and  500  mL/min,  respectively. 
The  power  reduction,  the  center  wavelength  redshift,  and  the  increase  in  FWHM  with  smaller 
heat  transfer  coefficients  are  evident. 
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Figure  61:  SPDs  of  the  LD  Array  measured  at  Heat  Transfer  Coefficients 

h  ,  =  56,250  W/(m2-K),  h2  =  80,500  W/(m2-K),  h3  =  112,650  W/(m2-K),  h4  =  119,900  W/(m2  K), 
and  hs  =  137,600  W/(m2  K)  with  If  of  80  A  and  Tiniet  of20°C 

The  measured  SPDs  of  the  LD  array  were  deconvoluted  by  the  proposed  method  to  investigate 
the  effect  of  the  heat  transfer  coefficients  on  the  SPD  of  individual  emitters.  The  results  of  four 
representative  emitters  (#1,  #4,  #8,  and  #12)  are  shown  in  Figure  62.  It  is  worth  noting  that  the 
maximum  power  of  the  edge  emitter  at  hi  was  the  lowest  among  all  edge  emitters,  but  the 
maximum  power  of  the  center  emitter  at  hi  became  the  highest  of  all  center  emitters.  This  is 
attributed  to  strong  coupling  between  the  junction  temperature  and  power  distributions. 
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Figure  62:  Predicted  SPDs  of  Emitters  #1,  #4,  #8,  and  #12  in  the  LD  Array  at  hi,  h2,  h3,  h4, 

and  hs 

In  order  to  put  this  into  perspective,  the  maximum  power  of  each  emitter  was  calculated  from 
Equation  (14).  The  results  are  plotted  for  hi  h2  and  hs  in  Figure  63  together  with  the  junction 
temperatures  predicted  in  Section  4.  The  current  competition  produces  significant  power 
variations  in  the  LD  array.  The  power  of  center  emitters  (from  Emitter  #8  to  #12)  is  almost 
constant,  but  it  decreases  rapidly  after  Emitter  #6. 

At  hi  =  56,250  W/(m2-K),  the  average  junction  temperatures  of  Emitter  #1  and  #12  are  47.4  and 
39.8  °C,  respectively,  and  the  maximum  power  of  Emitter  #1  is  about  28%  of  Emitter  #12. 
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At  hs  =  137,600  W/(m2-K),  the  average  junction  temperatures  of  Emitter  #1  and  #12  are  33.5  and 
29.9°C,  respectively,  and  the  maximum  power  of  Emitter  #1  is  about  54%  of  Emitter  #12. 

As  mentioned  earlier,  the  maximum  power  of  Emitter  #1  (edge  emitter)  at  hi  =  56,250  W/(m2  K) 
is  lower  than  that  at  hs  =  137,600  W/(m2  K).  However,  the  maximum  power  of  SPD  at 
hi  =  56,250  W/(m2-K)  increases  much  more  quickly  towards  center  emitters  compared  to  the 
case  of  hs  =  137,600  W/(m2-K).  As  a  result,  the  maximum  power  of  Emitter  #12  (center  emitter) 
at  hi  =  56,250  W/(m2-K)  is  higher  than  that  at  hs=  137,600  W/(m2-K). 


The  results  clearly  show  that  the  optical  power  ratio  between  center  and  edge  emitters  becomes 
larger  as  the  heat  transfer  coefficient  becomes  smaller.  It  is  known  that  the  power  efficiency 
reduction  of  each  emitter  caused  by  higher  junction  temperatures  shown  in  Figure  63  is  not 
significant. [61,  63,  64]  Therefore,  the  larger  ratio  at  the  lower  heat  transfer  coefficient  (i.e.,  the 
larger  temperature  variations)  is  mainly  attributed  to  the  current  competition.  The  larger  power 
ratio  between  emitters  are  expected  even  at  high  heat  transfer  coefficients  when  operating 
currents  much  higher  than  80  A  are  used. 


The  proposed  method  is  applicable  to  LD  arrays  with  higher  fill  factors  (up  to  90%)  if  the  SPD 
of  center  emitters  is  determined.  This  method  can  also  be  employed  to  deconvolute  the  SPD  of 
the  LD  array  at  different  operating  currents.  It  is  important  to  recall  that  the  SPD  profile  of  the 
center  emitter  are  altered  by  the  operating  currents  since  the  gain  spectrum  changes  with  the 
carrier  density  [47],  The  normalized  SPD  of  the  center  emitter  must  be  determined  at  a  given 
operating  current  for  successful  deconvolution. 
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Figure  63:  Junction  Temperature  and  Power  of  each  Emitter  in  the  Left  Half  of  the  LD 

Array  at  hi,  h2,  and  hs 

(If  =  80  A,  and  Tmet  =  20  °C) 
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3.3.5  Summary 


A  novel  method  was  proposed  to  predict  the  SPDs  of  individual  emitters  in  high  power  LD 
arrays.  The  objective  was  achieved  by  deconvoluting  the  SPD  of  an  LD  array  while  taking  into 
account  the  thermal  cross  talk  effect  as  well  as  the  current  competition  effect.  A  commercial 
water-cooled  LD  array  was  used  to  implement  the  proposed  method.  The  SPDs  of  individual 
emitters  in  the  LD  array  were  deconvoluted  successfully  at  different  cooling  conditions.  The 
results  indicated  very  strong  coupling  between  the  junction  temperature  and  power  distributions. 
The  comparison  between  the  predicted  SPD  and  the  experimentally  measured  SPD  showed  an 
excellent  agreement  in  both  shape  and  magnitude,  which  corroborated  the  validity  of  the 
proposed  method.  The  proposed  method  can  be  employed  to  improve  the  packaging  structure 
and/or  to  optimize  the  cooling  conditions  for  enhanced  pumping  efficiency  of  laser  diode 
pumped  solid-state  lasers  and  fiber  lasers. 

3.4  Lifetime  Prediction  of  LD/Microcooler  Subassembly  based  on  PoF  model 

In  this  section,  the  crack  propagation  in  the  die  attach  caused  by  the  thermal  fatigue  and  its  effect 
on  the  SPD  change  are  demonstrated  for  the  PoF-based  reliability  assessment. 

3.4.3  Crack  Propagation  Model  of  Die  Attach  for  LD  Bar 

The  viscoplastic  behavior  of  solders  should  be  used  to  predict  the  crack  propagation  in  the  die 
attach.  The  viscoplastic  behavior  of  solders  can  be  modeled  by  Anand’s  model  [65].  It  is  a 
unified  model  determined  directly  by  combining  both  rate-dependent  (creep)  and  rate- 
independent  inelastic  strains  (plasticity)  into  a  viscoplastic  strain  term.  Anand’s  model  requires 
inputs  for  the  nine  constants  and  can  be  expressed  as: 


: 


More  technical  details  of  the  Anand  constants  can  be  found  in  Ref.  [65],  Table  8  shows 
constants  of  Anand  model  for  SAC305. 
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Table  8.  Constant  of  Anand  Model  for  Indium  and  SAC305 


Parameter 

Name 

Sn3.0Ag0.5Cu  [68] 

A  (1/s) 

Pre-exponential  factor 

17.994 

Q/k(l/K) 

Activation  energy/Universal  gas  constant 

9970 

% 

Multiplier  of  stress 

0.35 

m 

Strain  rate  sensitivity  of  stress 

0.153 

hO  (MPa) 

Hardening/softening  constant 

1525.98 

a 

Strain  rate  sensitivity  of  hardening  or 
softening 

1.69 

SO  (MPa) 

Initial  value  of  the  initial  variable 

2.15 

S  (MPa) 

Coefficient  for  deformation  resistance 
value 

2.536 

n 

Strain  rate  sensitivity  of  saturation  value 

0.028 

Darveaux’s  approach  (energy-based  approach)  has  been  widely  accepted  to  predict  the  crack 
propagation  of  solders  caused  by  the  thermal  fatigue  [69].  The  energy-based  approach  utilizes  a 
finite  element  analysis  to  determine  the  inelastic  strain  energy  density  accumulated  per  each 
thermal  cycle.  The  strain  energy  density  and  the  crack  growth  data  are  used  to  predict  the 
number  of  cycles  to  initiate  and  propagate  the  cracks  through  a  solder  joint.  The  model  can  be 
expressed  as: 


da 


dN. 


£4 


where  No  is  the  crack  initial  life,  NP  is  the  crack  propagation  life,  a  is  the  length  of  crack.  The 
model  constants  of  Ki,  Ki,  K3,  and  K\  are  determined  empirically  from  the  test  data.  In  the  model, 
the  lifetime  becomes  the  sum  of  the  crack  initial  life  and  the  crack  propagation  life  (N  =  No  +  NP). 
A  Wave  is  the  average  inelastic  strain  energy  density,  which  can  be  calculated  using  Anand’ s 
model  [65].  The  averaged  inelastic  energy  density  change  per  thermal  cycling,  A  Wave,  is  defined 
as 


AW  - 


eiem 


TK 


ele. m 


eiem 


where  A Weiem  is  the  inelastic  strain  energy  density  change  per  thermal  cycling  of  each  element  in 
the  finite  element  model;  and  Veiem  is  the  volume  of  each  element. 
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In  the  model,  the  lifetime  becomes  the  sum  of  the  crack  initial  life  and  the  crack  propagation  life 
(TV/-  No  +  A/p).  The  life  time  can  be  expressed  as 


JL.-1 


a 


Table  9  shows  the  crack  growth  correlation  constants  of  the  SAC305  determined  based  on  the 
experimental  results  and  numerical  simulations,  respectively. 


Table  9.  Fatigue  Constants  for  SAC305  [70] 


Ki  (cycles/MPaK2) 

K2 

K3  (m/cycle/MPaK4) 

K4 

37.97 

-2.80 

1.4E-6 

1.16 

Figure  64  shows  the  FEM  model  of  the  half  model  and  the  die  attach.  The  die  attach  has 
thickness  of  5  pm  and  was  simulated  by  using  three  layers  of  elements. 


(a) 


Figure  64:  (a)  Half  Model  for  the  Solder  Model  and  (b)  Top  View  of  the  Die  Attach 
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Figure  65  shows  the  thermal  cycle  loading  conditions  used  in  the  study.  The  average  inelastic 
strain  energy  density  was  calculated  subjected  to  the  passive  thermal  cycles  with  the  military 
standard  condition  (-55°C  to  85°C,  ramp  rate  10°C/min,  dwell  time  15  minutes,  58  min/cycle) 
[66], 


Time  (min) 


Figure  65:  Thermal  Cycle  Loading  Conditions 


Three  different  crack  lengths  of  0.25,  0.75,  and  1.25  mm  were  considered.  Figure  66  shows  the 
three  crack  lengths  from  the  edge  of  the  die  attach.  The  crack  length  0.25  mm  does  not  reach  the 
first  edge  emitter.  The  crack  length  of  0.75  mm  penetrates  the  first  edge  emitter.  The  crack 
length  of  1.25  mm  penetrates  the  first  and  the  second  edge  emitters.  Figure  67  shows  the  plastic 
energy  density  distribution  of  the  die  attach  at  the  end  of  the  4th  cycle  used  for  fatigue  life 
calculation. 


0.25  mm  (5%  crack) 
0.75  mm  (15%  crack) 
1 .25  mm  (25%  crack) 


Figure  66:  Crack  Length  Information  from  the  Edge  of  the  Die  Attach 
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Figure  67:  Plastic  Energy  Density  Result  at  4th  Cycle 


Table  10  shows  the  number  of  cycles  required  for  the  three  crack  lengths.  The  number  of  cycles 
to  initiate  the  crack  (Ni)  was  3217.  The  total  number  of  cycles  (Nf)  for  the  SAC305  was  4340, 
6596,  and  8832  for  the  crack  length  of  0.25,  0.75,  and  1.25  mm,  respectively. 


Table  10.  Cycles  for  Crack  Length  with  SAC305 


Elements  for 
average 

Averaged  accumulated 
Strain  Energy  Density 
per  Cycle  AW  (MPa) 

Crack  length 
(mm) 

Ni 

(Cycles) 

NP 

(Cycles) 

Nf 

(Cycles) 

Whole  top 
layer 

0.205 

0.25 

3217 

1123 

4340 

0.75 

3217 

3379 

6596 

1.25 

3217 

5615 

8832 

3.4.2  Junction  Temperature  Distribution  Change  from  Crack  Propagation 


The  junction  temperature  distribution  was  predicted  for  4  different  cases:  the  crack  lengths  of  0 
(without  crack),  0.25,  0.75,  and  1.25  mm.  Figure  68  shows  the  loading  and  boundary  conditions 
used  in  the  thermal  model.  A  total  heat  power  of  30  W  was  applied  on  the  emitters  of  the  half 
model,  and  the  bottom  surface  temperature  of  an  LD  device  was  set  to  be  25  °C.  The  crack  area 
was  replaced  with  air  to  simulate  the  crack  propagation.  The  simulation  only  considered  the  heat 
conduction  due  to  the  negligible  heat  convection  and  radiation. 


w 


T0  =  25  °C 
at  bottom 


Figure  68:  Loading  and  Boundary  Conditions  for  the  Thermal  Model 
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Figure  69  shows  (a)  a  3-D  representative  view  of  the  temperature  distribution  of  the  half  model 
and  (b)  the  junction  temperature  distributions  of  emitters  without  the  crack  and  with  the  1.25  mm 
crack  in  the  die  attach.  The  junction  temperature  decreases  from  the  center  to  the  edge  emitter;  it 
also  does  from  the  emitting  side  to  the  back  end.  The  significant  junction  temperature  increase 
appears  at  the  edge  emitters  when  the  cracks  are  formed. 


Figure  69:  (a)  3-D  Representative  View  of  the  Temperature  Distribution  of  the  Half  Model 
and  (b)  the  Junction  Temperature  Distributions  of  Emitters  without  the  Crack  and  with 

the  1.25  mm  Crack  in  the  Die  Attach 


The  average  junction  temperature  of  each  emitter  without  a  crack  is  compared  with  those  with 
different  crack  lengths  in  Figure  70,  where  the  only  left  half  is  shown  due  to  the  symmetry.  The 
average  junction  temperatures  remain  virtually  unchanged  over  the  center  half  of  emitters  and 
then  drops  toward  the  edge  emitter.  As  expected,  the  junction  temperature  increased  significantly 
at  edge  emitters  with  the  cracks. 
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Figure  70:  Average  Junction  Temperature  of  each  Emitter  in  the  Half  of  the  LD  Array 
without  a  Crack  and  with  different  Crack  Lengths 

3.4.3  SPD  Change  from  Junction  Temperature  Distribution  Change 

The  SPD  change  was  predicted  for  4  different  cases:  the  crack  lengths  of  0  (without  crack),  0.25, 
0.75,  and  1.25  mm.  The  junction  temperature  distribution  obtained  from  the  thermal  model  was 
used  as  an  input  for  the  SPD  estimation. 

Three  Gaussian  functions  were  used  to  fit  the  normalized  SPD  of  the  center  emitter.  The  results 
are  shown  in  Figure  71  and  the  Gaussian  fitting  parameters  are  shown  in  Table  11. 


Wavelength  (nm) 

Figure  71:  SPD  Deconvolution  of  the  Normalized  SPD  of  the  Center  Emitter  using  Three 

Gaussian  Functions 
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Table  11.  Three  Gaussian  Functions  used  to  define  the  Normalized  SPD  of  the  Center 

Emitter 


Ai  (nm) 

A2  (nm) 

A3  (nm) 

807.980 

808.338 

808.608 

Wi  (nm) 

W2  (nm) 

W3  (nm) 

0.463 

0.285 

0.200 

A! 

a2 

a3 

0.655 

0.541 

0.549 

Figure  72  shows  the  normalized  SPDs  of  the  LD  array  measured  with  different  crack  lengths 
(0.25,  0.75,  and  1.25  mm).  The  wavelength  coefficient  was  0.3  nm/K  and  the  coefficient  of 
current  the  competition  was  5.9.  When  the  crack  size  is  less  than  0.25  mm,  the  total  SPD  of  LD 
bar  is  nearly  the  same  because  of  the  negligible  temperature  change.  It  should  be  noted  that  the 
edge  emitter  works  even  when  the  crack  is  propagated  because  the  emitters  are  connected  in 
parallel. 

The  center  emitter  of  the  SPD  also  showed  significant  shift  due  to  the  junction  temperature 
change  and  the  current  competition  effect.  When  the  crack  penetrates  to  one  emitter  (the  10th 
emitter),  the  emitter  temperature  increases  dramatically.  Most  of  the  power  will  be  loaded  at  the 
emitter  with  the  crack  due  to  the  current  competition.  Based  on  this  study,  it  is  expected  that 
there  is  no  gradual  degradation  due  to  the  thermal  fatigue,  but  the  fiber  laser  will  be  failed 
catastrophically  once  the  crack  penetrates  the  first  edge  emitter  due  to  the  significant  SPD  shift 
(coupling  efficiency). 


Figure  72:  Normalized  SPDs  of  the  LD  Array  Measured  with  different  Crack  Lengths 

(0.25,  0.75,  and  1.25  mm) 
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3.4.4  Summary 

The  crack  propagation  in  the  die  attach  due  to  the  thermal  fatigue  and  its  effect  on  the  spectral 
power  distribution  change  were  demonstrated  for  the  PoF-based  reliability  assessment. 
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LIST  OF  ABBREVIATIONS,  ACRONYMS,  AND  SYMBOLS 


ACRONYM 

DESCRIPTION 

AFRL 

Air  Force  Research  Laboratory 

COP 

coefficient  of  performance 

CTE 

coefficient  of  thermal  expansion 

DAQ 

data  acquisition  module 

FWHM 

full  width  at  half  maximum 

HTC 

heat  transfer  coefficient 

LD 

laser  diode 

LED 

light  emitting  diode 

MOSFET 

metal-oxide-semiconductor  field-effect  transistor 

PoF 

physics-of-failure 

SPD 

spectral  power  distribution 

TEC 

thermoelectric  cooler/cooling 

TIM 

thermal  interface  material 
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